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Chapter 1 
Introduction 



This Thesis deals with the one- dimensional (ID) Hubbard model [|1], 0, ^, which 
describes single-band interacting (correlated) electrons in a (ID) lattice. According 
to Lieb i: 

" The Hubbard model is to the problem of electron correlations as the Ising 
model is to the problem of spin-spin interactions; it is the simplest possible model 
displaying many "real world" features." 

The model depends on the parameter t, the hopping integral, characterizing the 
kinetic energy of the electrons when they hop between two adjacent lattice sites, 
and on the parameter U, the on-site Coulomb integral, characterizing the Coulomb 
interaction energy when two electrons occupy the same lattice site (Chapter H). 

Strictly one-dimensional systems are not found in real solid-state materials, how- 
ever there are quasi-one-dimensional materials in nature. Theoretically, a quasi-one- 
dimensional solid - the definition of which is quite broad 0- may be represented 
by a three-dimensional array of one-dimensional chains, such that the hopping in- 
tegrals between the chains are much smaller than the hopping integral along them 
1^. The correspondence with real materials - for which t and U can be estimated 
from infrared spectroscopy data R H, is then made by means of quantum chem- 



istry calculations [|T0| , where the hopping integrals and the Coulomb on-site energy 
are computed by means of the atomic (molecular) wave functions of the atoms 
(molecules) of the solid under consideration. 

Obviously, approximating real materials by a simple one- dimensional model is 
not enough in general to make a quantitative description of the solid properties. In 
spite of this, a description based on the parameters t and U contains many "real 
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world" features On the other hand, effects hke electron-phonon couphng 

and more general Coulomb interactions, for example, are needed to account for a 
more quantitative description of the physical properties of real quasi-one-dimensional 
materials [|Tl], |12|, |13|, [l^ . 

The systems described by the ID Hubbard model (and its generalizations) have 
narrow energy bands, and, therefore, the tight-binding approximation for indepen- 
dent electrons can be applied |]18 . 



1.1 Materials 

There are several classes of quasi-one- dimensional materials in nature which can 
can be described by ID electron models, the benzene molecule (CeHg) being one of 
the simplest systems to which the Hubbard model can be applied. [The molecular 
TT orbitals of the six carbon atoms can be thought of as a lattice with six sites 
populated by six electrons. The Coulomb interaction have been computed by Parr 
et al. [0 and are U = 16.93 eV, for the on-site Coulomb integral, and V = 9.027 
eV, for the next-nearest-neighbour Coulomb integral (Chapter 

In what follows we describe very briefly some of the most typical quasi-one- 
dimensional materials. 



1.) Inorganic Materials 



l.a) Potassium cyanoplatinates: a typical example is K2Pt(CN)4Bro.3.3H20 
which was synthesized during the 19th century. The chains are built by the super- 
position of the d orbitals of the Ft atoms. The ratio of the conductivity along the 
chains {a\\) over the conductivity perpendicular to the chains {(T±) is cr\\/cr± ~ 10^, 
and the electronic density n - the ratio of the number of electrons over the num- 
ber of lattice sites Na - is 0.85 |19|. The values of a\\ and a± are a measure for the 



atomic (or molecular) orbitals overlap along and perpendicular to the chain direc- 
tion, respectively [^. Therefore, cr\\/a± provides a measure of the relative values 
of the hopping integrals. 

l.b) Transition metalchalcogenides: MX3 and (MX4)„Y, with the metals 
M=Nb,Ta, the chalcogenides X=S,Se, and the halides Y=I,Br,Cl. The basic units 
of the chains in MX3 are MXg triangular prisms stacked on top of each other, and 
in (MX4)2Y are M2X4 octaedra. Typical examples are NbSea, TaSs, (TaSe4)2l, and 
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(NbSe4)2l, with a\\/a^ ~ (10 - 10^) and n ~ 1/4 [|T|]. 

l.c) Transition metal bronzes: Aq^sMoOs, with A=K,Rb,Tl. These materials 
are called blue bronzes due to their metallic blue brightness. The chains are formed 
by the MoOg octaedra, and, in these materials, o'\\/o-± ~ (10-10^) and n ~ 3/4 |TP[. 

l.d) Copper oxide compounds: typical examples are the ID antiferromag- 
netic Sr2Cu03 and SrCu02 materials |2l|. These systems can be modeled by the 



ID Hubbard Model with U » t and n = 1 (Chapter H). The Heisenberg ex- 
change interaction J = 4t^/f/ for these two compounds is J = (2200 ± 200) K and 
J = (2100 ± 200) K, respectively. The interaction between chains, J', obeys the 
ratio J/ J' ~ 10^ for SraCuOs and J/\J'\ ~ (5 to 10) for SrCuOa. Intimately related 
to the spin-chain compounds are the spin-ladder compounds. These latter systems 
can form two, three, or more leg ladders. Typical examples are SrCu203 (a two-leg 
spin 1/2 systems) and Sr2Cu205 (a three-leg spin 1/2 ladder) |^. An interesting 
result is that the even- leg spin 1/2 ladders present a gap in their excitation spectra, 
while the odd-leg spin 1/2 ladders do not. Upon doping (this means n = 1/2 — 6, 
where 6 is usually small when compared to 1/2) these ladder materials present both 
metallic and superconducting phases p2|, p3|. 



2.) Organic Linear Chain Compounds 

2. a) Polymers: the typical example is polyacetylene. It consists of a long chain 
of (CH)^; units with alternating single and double bonds between the carbon atoms 
fn\, O, p3. The values of the transfer and Coulomb integrals are taken to be 4t ~ 10 



eV and U ~ 6.5 eV, respectively, and n = 1 |11, 13, 24 



2.b) Stacks of organic molecules: these materials are made of planar or- 
ganic molecules stacked together. The organic molecules have vr molecular orbitals 
perpendicular to the plane of the molecule. These orbitals overlap forming a ID 
chain. A typical example is the incommensurate (the average number of carriers per 
molecule can not be expressed by a simple rational number) tetrathiafulvalena- tetra- 
cyanoquinodimethane - TTF-TCNQ - with At ~ 0.61 eV, U ~ 1.1 eV, and n ~ 0.55 
0, H. Other examples are the salts MsX, with M=TMTSF,TMTTF (where TMTSF 
denotes tetramethyltetraselenafulvaleno and TMTTF denotes tetramethyltetrathi- 
afulvaleno) and X=Br,PFe,C104. These salts are quarter filling systems (quarter- 
filled hole band: n = 1/4), with the t and U parameters being of the order t ~ 0.25 
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eV and f/ ~ 1 eV {th and tc, the interchain couphngs, are of the order of 25 meV 
and 1 meV, respectively) p5|. These materials can also present a superconducting 



phase [26 



There are still other realizations of quasi-one-dimensional systems like in het- 
erostructures (an example is the AlGaAl-GaAs-AlGaAl material) and in nanotubes 



(an example is the Bi_5C2+5N material) |^ |2^ . 

Many of the materials we have briefly described are not metals at zero tempera- 
ture. Usually they are metals at room temperature and undergo some type of phase 
transition when the temperature is lowered. Depending on the material and on the 
physical values of certain external parameters (like doping, external fields, pressure 
and cooling rate) they may, for example, undergo structural, metal-insulator, spin 



density wave, or superconducting phase transitions p6|, |19| 0, |3l|, |3^, |3^ . 

When one performs zero-temperature calculations and applies the results to 
finite-temperature data, no more than a quahtative agreement is to be expected. 
Nevertheless, the calculation of, for example, the zero-temperature optical conduc- 
tivity has for some quasi-one- dimensional materials been used with some success 
in describing the optical conductivity of these materials in the (finite temperature) 
metalhc phase 0, ^, [T^. However, due to its quasi-one- dimensional character, these 
materials can only be approximate by an one- dimensional model at moderate tem- 
peratures. At low temperatures 3D effects start to play a role. These effects may 
drive phase transitions in these materials, as those indicated above. 



1.2 Theoretical results 

As the Hubbard chain is the simplest possible model for the study of interacting 
electrons in a lattice, the theoretical study of the model is by it self of interest 
In this section we summarize very briefly some of the previous work on the ID 
Hubbard model. This includes its exact solution, as well as recent new results, some 
of which are presented in this Thesis. The ID Hubbard model is exactly solvable 
by the Bethe-ansatz technique, first introduced by Bethe in his solution of the 



ID spin 1/2 Heisenberg chain. Roughly speaking, the model has charge and spin 
excitations - usually called holon and spinon excitations - which are decoupled from 
each other. This feature is usualy taken as a signature of ID electronic correlated 
systems. The experimental observation of the holon and spinon spectra has recently 



been achieved by photoemission spectroscopy by Kim et al. [^] in the ID SrCu02 
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material (but not in the 2D related Sr2Cu02Cl2 cuprate). Moreover, these authors 
find that the Hubbard model in the limit U ^ At (the parameters used are t = 0.60 
eV and U = 7.2 eV) "agrees strikingly" with the data. 



Using the same method as Yang for a continuum multicomponent model |^ 
Lieb and Wu |^ computed the ground-state energy at half filling (n = 1) and zero 
magnetization (m = 0). They also found that for n = 1 the charge excitations have 
a gap - the Mott-Hubbard gap Amh^, for any value of the on-site Coulomb integral 
U. Thus the model is an insulator at half filling and zero temperature. Experimental 
evidence of ID Mott-Hubbard insulating behavior has been found in the quasi-one- 
dimensional HMTSF-TCNQF4, HMTTF-TCNQF4, and DBTSF-TCNQF4 organic 
compounds p|, |3^ (all these three compounds have n = 1, and HMTSF, HMTTF, 
and TCNQF4 denote hexamethylenetetraselenafulfalene, hexamethylenetetrathiful- 
falene, tetrafiuoro-tetracyano-p- quinodimethane, respectively). 

Shiba performed some numerical work for densities lower than half filling 
and computed the ground state energy and the magnetic susceptibility. In the limit 
U ^ 4t, it has been possible to compute analytically the ground-state energy for 
n < 1 The thermodynamics of the Hubbard chain has been described by 

Takahashi , who obtained the full solution of the Lieb and Wu equations . 



The spectra and the momenta of the low-lying excitations were determined 



by Ovchinnikov for the half-filled band case and by Coll for any filling |^ 



Ovchinnikov |^ also remarked that for U ^ At the spin spectrum is the same as 
for the anti-ferromagnetic Heisenberg chain |Q. This fact is an indication of the 
factorization of the Bethe-ansatz wave function into a free-spinless fermion Slater- 



determinant and into the anti- ferromagnetic Heisenberg wave function as dis- 
cussed later by Ogata and Shiba 

In the presence of a magnetic field, the charge and spin low-lying excitations 
transform into more exotic c and s excitations, and their energy spectra and mo- 
menta, as well as the low-energy thermodynamics and static properties have been 
analytically computed [46, 47, 4^ for electronic and magnetization densities such 
that n < 1 and m < n, respectively. The study of excitations of topological charac- 
ter in a finite-size chain has been made by Carmelo and Peres [^, In addition 
to the c and s modes, the ID Hubbard model has finite-energy charge (c, 7) and 
spin (5,7) excitations, which have been considered by Woynarovich [^l], These 
charge and spin excitations can be seen as bound states of c and s pseudoparticles. 
Their operator generators, excitation spectra, and momenta have been computed 
by Carmelo and Peres ||53[| . 
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Due to its involved form, the Bethe-ansatz wave function is not ideally suitable 
for the computation of operator mean values and correlation functions. These latter 
quantities can be obtained in some limits by combining the Bethe-ansatz solution 



with conformal-field theory |5|, 55|. These calculations require a simple representa- 
tion for the eigenstates of the model. In particular, the possibility of finding a set 
of operators (in second quantization) which diagonalize the Hamiltonian is a major 
step in that direction. An example of this general idea is the bosonization technique 



57| , where the g-ology Hamiltonians can be diagonalized by boson oper- 
ators. The low-energy correlation functions can then be computed by expressing 
both the fermionic fields and the eigenstates of these models in terms of the boson 
operators. 

A central problem in this model is the study of the transport properties: quan- 
tities like the charge and spin currents and the corresponding conductivity spectra 

The low-frequency properties of the ID Hub- 



are of importance [59, 30, pTl 32 



bard model have been studied both by means of bosonization techniques and 
by means of the pseudoparticle representation []63| , [6^ . However, the study of the 
finite-energy transport and spectral properties requires an operator representation 
for the Bethe-ansatz solution at all energy scales. 

I, Carmelo 



Following previous results for the low-energy Hilbert sub-space |65 



and Peres suceeded in writing the Hamiltonian, the momentum operator, and 
the generators for all 4^'' eigenstates of the Hubbard model in a suitable second 
quantized operational basis. This basis is defined by the set of operators ag^o,/3 
^g,a,7' which we call pseudohole and heavy-pseudoparticle operators, respectively 
(Chapter |^). In this basis the eigenstates of the ID Hubbard model have a simple 
Slater-determinant form. In its diagonalized form the model is written as an infinite 
sum of terms involving only forward scattering interactions among the pseudopar- 
ticles at all energy scales. 

The use of the pseudohole and heavy-pseudoparticle operational representation 
allows the evaluation of the charge and spin current mean values, as well as of the 
charge and spin transport masses of the pseudoparticle carriers [^, Moreover, 
combining this operational representation with a new generalized conformal-field 
theory ||6^, the optical conductivity of the ID Hubbard model around some finite 
energy values can be computed |7^. The set of these energy values, u = uq, corre- 
spond to the new conformal critical theories of {u — ujq) low energy. 

The present Thesis is organized as follows. In Chapter ^, a derivation of the 
Hubbard model is given and some of its well known properties are discussed. In 
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Chapter a detailed presentation of the operational representation of the model in 
terms of pseudoholes and heavy-pseudoparticles is given. The general expressions 
for the low-lying and finite-energy excitation spectra are provided in Chapter ^. 
In Chapter ^, the topological excitations in a finite size system are discussed. In 
Chapter ^ the charge and spin transport currents and the charge and spin transport 
masses are evaluated. Qualitative results on the optical conductivity of the Hubbard 
model are discussed and some Monte Carlo calculations are presented. Finally, in 
Chapter ^ we summarize and discuss future research directions. 
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Chapter 2 

The Hubbard Model 



2.1 The Hubbard model 

The theoretical study of electron correlations in sohds requires the use of a suitable 
Hamiltonian such that it includes, in some way, the Coulomb interaction among 
the relevant electrons of the sohd (those of the outer shells). If one is interested 
in both the thermodynamics and the transport properties, one must include the 
necessary external fields to the Hamiltonian. Response to an electric field requires 
the introduction of a vector potential A (taken to be small). In first quantization, 
the general Hamiltonian including Coulomb interactions reads 

1 

jEj-^. (2-1) 

where m is the electron mass, e is the electron charge (e < 0), and V{r) is the 
lattice potential, which will have the suitable symmetry (the effects of disorder will 
not be considered). The treatment of this Hamiltoninan is best accomplished in 
second quantization. We write the electron fields \E'^;(r) (cr and / are the spin 
and band labels, respectively) in terms of Wannier functions which perserve the 
lattice symmetry (for practical calculations, one uses atomic wave functions). These 
fields are localized at the lattice sites Ri and are of the form ^^^(r) = J2i4'{'^~ 

Ri)e ^"'^ '^"■''^^^^ '^^cl^i, where cj^^ creates one electron of spin a at the band I and 

— * J. 

lattice site Ri. The operators c]^ ^ and Cj^cr',i' obey the usual anticommuting fermionic 
relations and the exponential term in the field definition is required to mantain gauge 
invariance. As a first approximation we consider only one band and neglect all the 
band structure of the solid. We then omit the band index in what follows. In the 
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Wannier representation, the Hamiltonian (p.l[) reads 



i,j,l,ma,a' 



(2.2) 



where the hopping integral tjj- and the Coulomb integral Uij^i^m are given, respec- 
tively, by 



= - df(t>*{r- Ri 



2m 



+ V{r) 



(f-i?,). 



(2.3) 



and 



i,j,l,m 



dfdf'(j)*{f- Ri)(j)*{r' - Rj) 



Ri)<Pif-Rn 



(2.4) 



The hopping integral ( p. 3D gives the amplitude for an electron to jump from the 
atomic orbital centered at the lattice site Ri to the atomic orbital centered at the 
lattice site Rj. The Coulomb integral accounts for the electronic correlations among 
the electrons at the different atomic orbitals. Hubbard estimated [|l| that the most 
important contribution to the Coulomb integral comes from the term where two 
electrons are at the same atomic orbital, that is f/i,j,j,j = U (the order of magnitude 
of the bare U is roughly of 10-20 eV). The second most important contribuition 
to the interaction energy, Uij^ij, comes from the interaction of two electrons at 
nearest neighbor sites (the order of magnitude is roughly 2 to 3 eV). If the electrons 
belong to narrow free-electron bands, the most important contribution to the kinetic 
energy comes from hopping between nearest neighbour sites; this is the usual tight- 
binding approach. In view of this, the simplest Hamiltonian we can think of for 
describing correlated electrons in a regular solid is the Hubbard model. It consists in 
considering, out of all the hopping and Coulomb integrals, only the nearest neighbour 
hopping t and the on-site Coulomb repulsion U. The Hubbard Hamiltonian then 
reads ig § 



H 



Na 

-t E 

<i,j>,a 



^ ch 



(2.5) 



where < i,j > means summing over nearest neighbor sites i and j only, Na is 
the number of lattice sites, and fii^^ = c-^^Ci^a is the number operator at site i for 
a electrons. For the spin label, a, we have a =t, i when used has an index and 
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(J = ±1 otherwise. This Hamiltonian has only two parameters, t and U , and thus 
the only energy scale is U / 1. This ratio is present in all physical quantities. In what 
follows we consider c^^+i^ = Cio-, units such that h = 1, the lattice spacing a = 1, 
and the electron charge e = — 1. 

2.2 The SO(4) symmetry of the Hubbard model 

The SO (4) symmetry of the Hubbard model can be traced back to its spin- up/spin- 
down and particle-hole symmetries ^, 0. Let us define the number and magne- 
tization operators as iV = J2i,afi'i,o- and = — | Si,o- cr^i.o-! respectively. {N^ is the 
number of a electrons and = N^r + N^^-) If we perform the spin- up/spin-down 
transformation, cfg. c] the local number operator for a electrons transforms as 
'fi'i-a- This implies that the number operator, N, and the Hamiltonian ( p.5| ) 
remain unchanged, but that the magnetization operator transforms as ^Sf.. 
That is, under this transformation we go from the sector where N^r > N^a- to the 
sector where Nf, < A^_o-. This can be accomplished by the ladder operators 

^- = EsV.i, ^+ = EsV.T- (2-6) 

j 3 

On the other hand, if we perform the particle-hole transformation Ci^^ 
the local number a operator transforms as ^j o- (1 — ^j.-o-)- This transformation 
implies that the magnetization operator and the kinetic energy {A = 0) of the Hamil- 
tonian (|2.5|) remain unchanged, whereas the number operator and the interaction 
term of the Hamiltonian (|2.5| ) transform as 

N-^2Na-N, (2.7) 

and 

U E ^ ^ E + " UN , (2.8) 

i i 

respectively. Therefore, the interaction term can be cast in an invariant form if we 
write it as f/I]j(^ii.| — 1/2) (nj j — 1/2). This transformation then takes us from the 
sector where N < Na to the sector where N > Na- At half- filling the number opera- 
tor remains invariant. This change in the electronic density sector without changing 
the magnetization density can be accomplished by means of suitable operators. 
The simplest way for generating these operators is to perform the transformation 
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Ci^^ — > (— l)*cj| in the St, 5*^, and SI operators. This leads to the following set of 
operators 

S: = -llNa - E^-] ' Sl= E(-iy9TC,i , SI = E(-1)^'sV]t • (2-9) 

Both the set of generators {St, S^, S^} and {St, S^, S^} form two independent 
SU(2) algebras. The first set takes into account the spin part and the second set 
the charge part. If we write the Hamiltonian ( ^.51 ) as 

Hso(4) = -tJ2 clc^,. + Uj2if^,,^ - l/2)(^u - 1/2) , (2.10) 

the Hubbard model becomes SO (4) invariant 0, ^, 0. Note that the SO [4) sym- 
metry is not fully equivalent to the SU{2) ® SU{2) symmetry 0, 0. All the 
eigenstates |0) of the Hubbard model can be choosen as simultaneous eigenstates of 
H, S'^, S^, {S'^y, {S'^y, and the lattice translation operator. The commutators 

[Sl,N] = -2Sl, [Sl,S:] = -S^, a = c,s, (2.11) 

and 

[S%P] = -7rSl, {P = ^kcl^c,,^) (2.12) 

where P is the momentum operator, imply that if a state has electrons, 
magnetization {N^ — N^)/2, and momentum P, then the state 5*^10) has + 2 
electrons and momentum P + vr, and the state 5*^10) has magnetization {Ni — 
A^l)/2 + 1 and momentum P. These results introduce a enormous simplification in 
the description of the Hilbert space in any dimension. 

The inclusion of terms that break the SO (4) symmetry in a trivial way, like the 
chemical-potential term and the Zeeman-coupling term, do not change the simplicity 
of the SO (4) description. The inclusion in the Hamiltonian (|2.10|) of such terms 
reduces, in general, the symmetry of the model. The resulting Hamiltonian reads 

ii = HsoiA) + 2fif], + 2fioHS, . (2.13) 

For finite values of both the magnetic field and the chemical potential the symme- 
try of the quantum problem is reduced to U{1) Cg> f/(l), with 5*^2 and S'^z commuting 
with H. The values of fi and H determine the different symmetries of the Hamilto- 
nian (|2.13|) . When /i 7^ and H the symmetry is U{1) U{1), for fi = and 



2.3 Basic features of the 



19 



^ it is SU{2) O f/(l), when /i and = it is f/(l) ® Sf/(2), and at /i = 
and H = the Hamiltonian symmetry is SO (4). 

There are four f/(l) ® f^(l) sectors of parameter space corresponding to S'^ < 
and SI < 0, < and > 0, S^, > and < 0, and > and SI > 0. We 
follow the works [§, ^ and call these sectors (Ich), where 




(2.14) 



The sectors (—1, —1), (—1, +1), (+1, —1), and (+1, +1) refer to electronic densities 
and spin densities < n < 1 and 0<m<n, 0<n<l and —n <m<0, l<ri<2 
and < m < (2 — n), and 1 < n < 2 and —(2 — n) < m < 0, respectively. 

There are two (Ig) sectors of SU{2) ® U{1) Hamiltonian symmetry [and two (Ic) 
sectors of U{1) (E) SU{2) Hamiltonian symmetry] which correspond to S*^ < and 
5*1 > for /<, = — 1 and Is = +1, respectively, (and to 5*^ < and S'f > for Zc = — 1 
and Ic = +1, respectively). There is one S0(4) sector of parameter space [which is 
constituted only by the S'^ = (and /i = 0) and S^ = canonical ensemble]. 



2.3 Basic features of the Hubbard model 



The Hubbard model ( |2.5| ) [or its SO (4) invariant version ( |2.1CI| )] depends on two 



parameters, namely the hopping integral t and the Coulomb integral U. The model 
shows very different physical behavior depending on the relative strength of these 
two parameters and on the band filling n. In the case of ?7 = 0, the Hubbard model 
reduces to the simple tight-binding model 

Ho = -t Y: e-'^'^-^^-^M,.c,,. . (2.15) 

<i,j>,cr 



and a simple exact solution follows from |]I0| the introduction of momentum-space 
operators ct and c^^. These operators are related to the Wannier representation 

by cJcr = A^""^/^ X]j e^'^'^'ct^. The form of the band then depends on the type of 
lattice considered. For the one-dimensional case, with one atom per unit cell, the 
diagonalized Hamiltonian is 

^0 = -2t cos (k - -) cl^^Cfc,. . (2.16) 



c 



It is a simple matter to compute both the exact partition function and the conduc- 
tivity for the Hamiltonian ( p. 161) . This model is a metal for all electronic densities. 
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exception made to N/Na = 2, where it has two electrons per lattice site, and for 
N/Na = 1, when the system is fully polarized. 

Another limit where the model ( p.5|) has an exact solution is the atomic limit 
In this limit the hopping integral is zero and the electrons are localized at the 
lattice sites. Since the electrons have zero kinetic energy, the model represents an 
insulator for all values of U, possesing only two (flat) bands, one located at the 
chemical potencial /i and another at the energy fi + U. This can be put in a more 
formal way by means of the Green's function formalism. A very simple calculation 
of the one-particle Green's function at finite temperature (suitable to deal with the 
massive degeneracy of the model eigenstates in this limit), leads to 



'2.17) 

iu + fi iu + fi — U \ ' ^ 

where iu is the usual finite-temperature Matsubara frequency and {n^) is the elec- 
tronic occupation number at finite temperature for the single site Hamiltonian, 

The question now is whether the system undergoes or not a phase transition 
when the parameters varie from the tigth-binding U/t = limit to the atomic 
U/t —>■ oo limit, for some finite value of U/t. In the strong-coupling regime U ^ t, 
the states with two electrons in the same site have a much larger energy - of order U 
- relatively to those states with single-site occupancy only. Thus, for the half-filled 
band we expect to have an insulator in the large U limit, and we show below, within 
perturbation theory, that this is indeed the case. 

Let us then divide the Hilbert space of the Hubbard model into three sub- 
spaces: the subspace of single-occupied-site states {15")}, the subspace of one double- 
occupied site states {|Di)}, and the subspace of more than one double-occupied-site 
states Since U ^ t, the kinetic term of the Hamiltonian (|2.5| ) can be treated 

as a perturbation. Note that due to the massive degeneracy in each of the three 
considered subspaces, degenerate perturbation theory must be applied. Moreover, 
in order to extract information on the nature of the ground state of the model, we 
compute the average energy relative to the states belonging the subspace 

Let us consider the system away from half filling. All states belonging to the 
subspace {15")} have zero Coulomb energy. Then, up to second order in perturbation 
theory, the energy of a given state belonging to the subspace is given by 

E\s) = {S\f\S) - ^ Y.iS\T\D,){D,\f\S) , (2.18) 

^ \Di) 
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where T = —tJ2<i,j>,cT cl_crCj,a- The states belonging to the subspace {\D„)} have not 
been considered because the operator T does not connect that subspace with {15)}. 
The operator P^,^ = J2\Di) \Di){Di\ can be regarded as the projection operator for 
the subspace Thus, up to this order in perturbation theory, the operator 

Pdi can be replaced in Eq. ( p.l8|) by I]j rii,!^^ j. It should be emphasized that this 
substitution is valid only if these calculations are restricted to the subspaces {\S)} 
and {|Z?i)}. From Eq. ( |2.18| ), one then arrives to an effective Hamiltonian valid for 
the subspace {\S)} that reads 



Tjeff. _ p 



a, a' l^i,m 



(2.19) 



where Pg = J2\s) \^){^\ is the projector onto the subspace The third term in 

the rhs of Eq. ( p.l9|) is usually neglected and the remaining two terms are called 
the t — J model 0]. At half filling, the first and the third terms in the rhs of Eq. 
( p.l9| ) vanish and we are left with another effective Hamiltonian 



m = ^T.\S^-S,-'^]. (2.20) 

with J = At'^/U. This is the Heisenberg anti-ferromagnetic model. The second term 
in the rhs of Eq. (|2.19|) was rewritten by using the electronic representation for the 
spin operators defined by S^, S^, and St- 

The above results show (assuming that the perturbative treatment holds) that for 
strong electronic correlations {U ^ t) the system is an antiferromagnetic insulator at 
half filling. Due to the restrictions imposed on electronic movements by the reduced 
dimensionality, for one-spatial dimension the Hubbard model is an insulator at half- 
filling for any positive value of U ||12[. In one dimension the system is not a simple 
Fermi liquid. The most difficult problem is the study of the intermediate-?/ case and 
of how the model properties change as the ratio U/t goes from small to large values. 
The study of this regime requires non-perturbative methods. In the following we 
consider the case of the ID lattice where the exact Bethe-ansatz solution is available. 
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2.4 Basic ideas on the Bethe ansatz 

From here on, we will be dealing with the one-dimensional Hubbard model. Follow- 



ing Yang's solution of the one-dimensional many-body problem with a delta- 



function interaction, Lieb and Wu ||T^ computed the ground state of the Hubbard 



model. Later, Takahashi |]I4| gave the complete solution of the Lieb and Wu Bethe- 
ansatz equations. The basic ingredients of the Bethe ansatz are: the simple topology 
of the ID space, the continuity condition at the same point of the lattice, the bound- 
ary conditions, and Pauli's principle. However, the implementation of these simple 
concepts in order to achieve the solution of the general eigenvalue problem is a 



task of considerable complexity In this section, the Bethe ansatz technique is 
illustrated by the solution of the two-electron problem in a chain of Na sites. 

The eigenstates |\E') of the two-electron problem can be written as a superposition 
of Wannier states cl^^a2^xi,ai\^) 

1^) = E E ^(a;i,ai;a;2,a2)4,.,4,.jO), (2.21) 

where |0) is the electronic vacuum and Xj stands for the position of the electron i. 
Obviously, the problem consists in determining the coefficients ip{xi,ai;x2,cr2). The 
starting point is the eigenvalue equation Hl"^) = El^^) which imposes the following 
relation among the coefficients ipixi, ai, X2, o"2) 

- t['ip{xi + 1, (Ti; X2, (T2) + ipixi - 1, (Ji; X2, (T2) + 

+ ll){Xi, Oi]X2 + l-,CT2) + ^{xi, ai;x2 + 1, (T2)] + 

+ U6xi,X26ai,-a2^{xi,ai;x2,a2) = E'ip{xi,ai;x2,(J2) ■ (2.22) 

The spatial dimensionality of the system allows the definition of two regions, I and II, 
for which xi < X2 and xi > X2, repectively. The electronic antisymmetric character 
under the interchange of positions implies that 

^'/(xi, di; X2, da) = -i'ii{x2, ^2; Xi, cxi) . (2.23) 

The next step is to choose a simple representation for the anti-symmetric property 
([^.23|). We introduce the Bethe ansatz 



M^u X2, a2) = A(pi, ai; p2, a2)e*(^^^^+^^^^) - A{p2, a,;p,, a2)e'(p^^^+P^-^^ (2.24) 



2.4 Basic ideas on 
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and 



It is clear that the interchange (a;i, ai) ^ (,T2, o"2) in i^ii^i, ai, X2(T2) or ipn{x2cr2, Xi, cxi) 
agrees with the antisymmetry condition ( |2.23| ). If xi ^ X2, introduction of both 
ipi{xi,ai;x2<J2) and 4'ii{x2<J2,Xi,ai) in Eq. ( p.22| ) leads to = — 2tcos(pi) — 
2icos(p2) (we are considering the lattice spacing a equal to 1). 
At the same lattice point, xi = X2 = x, we must have 

^i{x, CTi; X, 0-2) = i^n{x, (T2] x, (Xi) , (2.26) 

which is the continuity equation. Combining Eq. ( p.26|) with the eigenvalue equation 
for the two electrons at the same site and imposing the same functional form, E = 
—2tcos{pi) — 2tcos(p2), for the energy eigenvalue, after some algebra the following 
relation among three of the amplitudes A{pj , Pm, CTn) (with j,i,m,n = 1,2) is 
obtained 



/ A{p2, ai;pi, ^2) \ _ 1 / ic ki2 \ / A{pi, ai;p2, ^2) \ ^7) 
\A{p2,a2;pi,ai) J ki2 + ic\ki2 ic ) \A{pi,a2;p2,cri) J 

Here ki2 = sin(pi) — sin(p2) and c = U/ (2t). Solution of the problem in closed form 
involves the imposition of periodic boundary conditions (other boundary conditions 
are also possible). Since we have two electrons, we must impose periodic boundary 
conditions on both xi and X2 (some care is needed in taking into account the regions 
I and II). These are 



and 



Region I : xi < X2 ^ Xi + Na > X2, 

ipiixi, CTi; X20-2) = ^//(xi + Na, cTi; X2, 0-2) 



Region II : Xi > X2 ^ X2 + Na > Xi, 

^ii{xi, ai; X2a2) = ^lixi, ai;x2 + Na, as) 



(2.28) 



(2.29) 



Combinig Eqs. ( p.24| ) and ( |2.25|) with Eqs. ( p.28| ) and ( p.29| ), we arrive to the two 
following relations 
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\A{p2,a2;pi,ai) J \A{pi,ai;p2,cr2) J 

and 

^iP2Na f A{pi,cri;p2,cr2)\ ^ / A{p2,(T2;pi,ai)\ 
\A{pi,cr2;p2,cri) J \A{p2,ai;pi,a2) J ' 
By combining Eqs. ( [^.30| ) and ( |2.27| ), we end up with the simple eigenvalue problem 

e-'P'^-\A) = R\A) , (2.32) 

and equivalently for the eigenvalue e~*^^^°. The state \A) and the operator R are 
written as 



1^) 



' A{pi,ai]p2,(r2) 
^A{pi,a2]P2,cri) 



R 



1 (0 1' 
ki2 + ic I 1 



ic ki2 
ki2 ic 



(2.33) 



The solution to this simple problem can be cast in the form 
gip.TV, _ sm{pj) - A + ic/2 

and 



sin(pj) — A — ic/2 



J 



1,2, 



(2.34) 



sin(j9j) -A + ic/2 _ ^ 



(2.35) 



fj^ sm{pi) - A - ic/2 

where A = sin(pi) + sin(j92)- These are the Bethe-ansatz equations for one up-spin 
and one down-spin electron. 

Finally, by considering a magnetic flux through the ring, the general Bethe- 



ansatz equations become |T2|, |T6 

N 



n 



I sin(A;j) - Ks + iU/A 



\ sm{kj) - As- iU/A 



{j = l,...,N) 



(2.36) 



and 



/i, sin(fc,) - A^ - tU/A ~ A/3 - A5 - iU/2 



{6 = 1,...,N^). (2.37) 
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Chapter 3 



Algebraic Representation for the 
Hilbert Space of the ID Hubbard 
Model 

3.1 Introduction 

In this chapter we present the algebraic solution for the ID Hubbard model, which 
shows some basic similarities with the algebraic solution of the isotropic harmonic 
oscillator. In order to clarify the basic ideas underlying the problem, we begin by 
presenting in the following the main results for the algebraic solution of the harmonic 
oscillator. Afterwards, we discuss the general similarities and differences between 
the two cases. 

The isotropic harmonic oscillator has SO (3) symmetry. The SO (3) group can be 
defined by the following algebra envolving the angular momentum operators 

[L+,L,] = -L+, [L_,L,]=L_, [L+,L_] = 2L,, (3.1) 

where L = {L^, Ly, L^) and L± = ±iLy. Moreover, [L'^,L] = and [H,L] = 0, 
where H represents the isotropic harmonic oscillator Hamiltonian. The commutators 
( p.l| ) imply that L^, L^, and H can be simultaneously diagonalized. In second 
quantization, H, L^, and can be written as {h = 1) 

H = — [3 + alao + a\ai + a^ia_i ) , (3.2) 
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Lz = -{a\ai — a^_ia_i) , (3.3) 

2 2 2 2 2 

= L_L+ + LI + L,, (3.4) 

with 

L„ = itti , L,= a\a_i . (3.5) 
-52 ^22 

The operators Oq, a\, and a^i are bosonic operators and the angular momentum op- 

2 2 

erators are written in the Schwinger representation The eigenstates |ni , no, n_i) 
of H, Lz, and L^, corresponding to a given n = ni + no + n_i, a given i = n — 2no, 
and given m (the eigenvalue of L^), read 



2 2 



^{n-2no+m)^^t_J2(n-2no)(^t)4no 

\ni,no,n_i) = ' |0,0,0)- (3.6) 

y(4no)!y(2n - 4no)!v/(n - 2no + m)! 

For a given value of n, the quantum number no assumes the values 0,1,2, ... ,n/2, 
for n even, and 0, 1, 2, . . . , (n — l)/2, for n odd [the number no is used to generate, 
for example, all the LWS (m = —£, see below) for a given n]. The states of the form 
|ni , riQ, 0) are called highest-weight states (HWS's), since L+|rai , uq, 0) = 0, and the 
states |0,no,n_i) are called lowest-weight states (LWS's), since L^\0,nQ,n_i) = 0. 
The states such that L±\ni,no,n_i) ^ are called non-LWS's or non-HWS's. To 
go from the LWS of a given tower, characterized by the angular momentum number 
i = n — 2no, to the respective HWS (or vice- versa), one applies the rising operator 

(or the lowering operator L_). For a given i, the number m can have the values 
m = —i, — I + 1, ... ,1 — 1,1. All the states in the £ tower differ only in the 
eigenvalue, and all the states with a given n have degeneracy (n^ + 3n + 2)/2. 

From the above analysis, we see that the isotropic harmonic oscillator is char- 
acterized by three types of particles: the aj bosons, which do not contribute to the 

angular momentum operators, and the a\ and a^_x bosons. The latter two parti- 

2 2 

cles describe the towers of angular momentum. Note that the HWS's (LWS's) have 

only at (a^i) bosons. The non-LWS's or non-HWS's have both types of particles. 

2 2 

The ladder operators L± act upon the ±| bosons interchanging the ±| quantum 
numbers. 

A similar algebra occurs in the ID Hubbard model. Obviously, the algebraic 
description of this model is much more involved, for we are dealing with a many- 
body system with discrete translational symmetry. The Bethe-ansatz solution of 
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the Hubbard model refers usually to the symmetry sector (—1, —1) [Q, ^ Q (Sec. 
2]2|). In this sector it has been proved by Essler et al. 0] that the Bethe- 



ansatz eigenstates are LWS's of the S0(4) algebra, that is S'^\il)) = 0. We can 
then construct the towers of spin and eta-spin by acting onto the states {ip) with 
the ladder operators and reach the respective HWS's. This procedure generates 
all the 4^° eigenstates of the Hubbard chain, which span the whole Hilbert space 
associated with all the nine symmetry sectors of the model. 

In the algebraic solution of the Hubbard chain the eigenstates of the model are 
described by the set of operators a^ a^p and fej^a^^ which refer to the quantum objects 
that we call pseudoholes and heavy-pseudoparticles, respectively ||, ||, |10|. These 



operators obey the usual anticommutation relations. The quantum numbers that 
label the operators ctq^Q,,/? ^'^d ^J^Q^-y are the following: a) the q is the pseudomo- 
mentum and is to be chosen out of the available pseudo-Brillouin zone values (to 
be defined later); b) there are two colors of pseudohole and heavy-pseudoparticle, 
which are called a = c and a = s; c) these different types of heavy pseudoparticles 
can populate different 7 bands (7 = 1,2,..., 00). In the case of the a, P pseudoholes, 
the band index number is 7 = (this was omitted for convenience of notation); d) 
the /3 = ±1/2 is related to the S0(4) algebra. 

In the ID Hubbard model the Bethe-ansatz eigenstates can be classified as LWS's 
I and LWS's II [U ||, 0]. The former have neither heavy-pseudoparticle nor = 
1/2 pseudohole occupancy and all the pseudoholes which populate the a,0 
band have (3 = —1/2. On the other hand, the latter eigenstates have both heavy- 
pseudoparticles and P = ±1/2 pseudoholes. The SO (4) generators involve only 
the Og^Q^/j pseudohole operators, as the SO (3) generators, in the isotropic oscillator, 

involve only the 1 boson operators. The two colors a = c and a = s are directly 

^2 

related to the charge and spin degrees of freedom, respectively. The generators of 
the SO (4) algebra are given, in the generahzed Schwinger representation, by H, |T0| 



= -y^(a^ ia_„i— ia„^ 1), S'? = y^a^ _i_iO„„-ri. (3.7) 

^ q q 

The spin and eta-spin "angular" quantum numbers S*" are given by 

-1 00 

5" = ^(iV„'-E27iV«,7), (3-8) 
^ 7=1 

where A^"^ = J^fs^a^is: ^0,13 is the number of a,P pseudoholes, and Na,^ is the 
number of heavy a, 7 pseudoparticles. As in the S0(3) algebra, the ladder operators 
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interchange only the (3 = ±1/2 quantum number. Thus, as in the case of the 
harmonic oscillator, all the pseudoholes have (3 = —1/2 in the LWS's I and (3 = +1/2 
in the HWS's I. The non-LWS's I (and non-HWS's I) have both p = +1/2 and 
P = —1/2 pseudohole occupancy. 

The number of discrete pseudomomentum values in each a, 7 band (a,7-band 
Fock-space dimension) is given by the Bethe-ansatz solution M and reads 



4,0 = Na (3.9) 

00 

4,7 = A^a+A^a,7- Eb + 7'-|7-7'l]A^a,7', «, 7 7^ C, . (3.10) 

7'=1 

In the summations and products over q the pseudomomentum q takes on values 



271 

where in contrast to the usual momentum, J"'^ are consecutive integers or half-odd 
integers for N^^^ odd and even, respectively 0, ^ . Here 

A^c,o = E ^ - E Nc,, = ^ - iVs,o - E ^-,7 , (3.12) 

OL 7=1 0,7=1 

and 



for values of a, 7 other than c, 0. It follows that for each a, 7 band, q^^ < q < q^^, 
with the limits of the pseudo-Brillouin zones given by 

qif = ±n[l-^], (3.14) 

for Ncfi even and 

gS^ = vr, g(-) = -7r[l--|-], (3.15) 

or 

li'S =^{l-^], gS = -vr, (3.16) 
for Ncfi odd and simply given by 
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qi±) 



TT 



(3.17) 



for all the remaining a, 7 bands. 

In the LWS's II (or HWS's II) the situation is somewhat more complex, and 
has no analogon in the case of the harmonic oscillator. The origin of these states 
is in the lack of continuous translation symmetry. In these states, which can be 
interpreted as bound states of a,0 pseudoparticles, we have both /3 = 1/2 and 
/? = —1/2 pseudohole occupancy. The general second quantized representation for 
all Hubbard-chain eigenstates is 



|V';{<^},{iV„,,}) 



7=1 



q,a, 



1^) 



(3.18) 



where the symbols {N^^jj} and {Na^-y} are abbreviations for the sets {N^^i, Nl^_i, 
N^i,N^_i} and {iV,,i, iV,,oo, iV.,!, ^s,-}, respectively, C = UaiN^Wl^M 

"'2 '2 "'2 

= n/3 0(^a,/3 ~ Z^^^i 7^0,7)5 and the ladder operators S*" are defined in Eq. 
(|]7D. Note that the states ( glsD have a very simple form and present a remarkable 
similarity with the states ( p.6|) . 

As in the harmonic oscillator, the Hubbard model Hamiltonian (and other op- 
erators) can be written in terms of the operators and feg^Q^-y and the energy 
of the LWS's and HWS's is independent of the P quantum number. In general, 
the non-LWS's (non-HWS's) and the LWS's II (HWS's II) have an energy gap rel- 
atively to the corresponding ground state (which is always a LWS's I or HWS's I). 
In the remaining sections, we present the details of our algebraic solution for the ID 
Hubbard model. 



3.2 Pseudoholes: the (3 quantum number in LWS's 
I/HWS's I 

As was referred above, the description of the non-LWS's, not included in the Bethe- 
ansatz solution requires the introduction of the (3 = ±1/2 pseudohole numbers. We 
start the analysis with the study of the LWS's I. The Bethe-ansatz solution in the 
(-1,-1) sector (Sec. |2^ ) gives for the numbers of a,0 pseudoparticles, Nafi, and for 
the total numbers of permitted orbitals, d^fi (number of states in the respective 
pseudo-Brillouin zone), the following values P, ^ 
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iV,,o = iV, Ns,o = \{N,,o-N^ + N^) , (3.19) 

4,0 = Na , 4,0 = N,^o - • (3.20) 

The use of the S'0(4) symmetry (Sec. p.2| ) allows one to obtain these numbers in 
the remaining three U(1)®U(1) symmetry sectors. The number of pseudoholes in 
the a,0 band is given by A^^ = 4o "~ ^ao- The corresponding results for the 



LWS's I/HWS's I in the four U(1)®U(1) symmetry sectors are listed in Table [O . 
The numbers and in Table ^A] are those of the respective (/c, h) sector, and 
are related to those in different sectors by the particle/hole and spin-up/spin-down 



symmetries discussed in Sec. 2.2 



Defining A^^e as the electron number in a given eta-spin LWS, the respective 
eigenvalue of 5*^ is = —S'^ = —^{Na — N^,) in terms of electrons, or S*^ = — S''^ = 



— ^A^ _i (see Table 3.1), in terms of pseudoholes. On the other hand, the corre- 

^ C, 2 

spending HWS has 2A'a — A'e electrons and = S'^ = |(A^a — Ne), or, equivalently, 
= S'^ = IA^^i (see Table |3Tl|) . Obviously, equivalent results hold for the spin tow- 



ers. It follows from the results of Table ( |3.1| ) that the number of pseudoholes in any 





(-1,-1) 


(-1,1) 


(1,-1) 


(1,1) 




N 


N 


2Na-N 


2Na-N 








Na 


Na 




Na-N 


Na-N 








AT'^i 








N-Na 


N-Na 




Ni 




Na-N^ 


Na-N^ 


4,0 






Na-N^ 


Na-N^ 


* -2 







N^-Ni 





AT^i 





Ni-N^ 





Ni-N^ 



Table 3.1: Values for the total number of pseudoparticles Nafi, total number of 
accessible orbitals da o, and total number of pseudoholes A^'* i for LWS's I/HWS's 

' '^"^ -"-2 

I in the four sectors of symmetry U(1)®U(1). The numbers A^ and N^j are those of 
the respective (/c, h) sector. 



(/c, Is) sector is equal to 215*" | — the eigenvalues of S" — , defined in Sec. in terms 
of electronic operators and in Eq. ( |3.7| ) in terms of pseudohole operators. Then, in 
going from a LWS, belonging to a given sector, to the respective Sa HWS, we 
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have to apply the ladder operators 5*" a A*"^ number of times, with A^^ = Na — 
for a eta-spin tower and A^^ = A^^i — A^^ for a spin tower (the numbers and A^j^ 
are those for the respective LWS). Hence, as in the isotropic harmonic oscillator, 
we have in the LWS's I occupancies of /? = —1/2 pseudoholes only. On the other 
hand, the HWS's I show only /3 = 1/2 pseudoholes occupancy. In between we have 
non-LWS's or non-HWS's whose occupancies refer to mixtures of both j3 = —1/2 
and P = +1/2 pseudoholes. The number of pseudoholes is such that it gives the 
correct values for the eigenvalues and the correct number of states for a given 
S"' tower. 



3.3 Heavy-pseudoparticles: the p quantum num- 
ber in LWS's II/HWS's II 

In the LWS's II/HWS's II there is both pseudohole and heavy-pseudoparticle occu- 
pancy. In the (—1,-1) sector, the Bethe-ansatz solution gives the following results 



for the numbers of a, pseudoparticles and accessible da,o orbitals P, |T0| 



oo oo oo 

A^,,o = iV - 2 ^ 7Ar,,^ , Ar,,o = iV^ - E - E (7 + 1) > (3-21) 

7=1 7=1 7=1 

and 

oo 

4,0 = Na , 4,0 = N-2N^ + Ns^o + 2 E 7A^.,7 , (3-22) 

7=1 

respectively. From the above two equations, it follows that the number of a, 
pseudoholes populating the LWS's II is given by 



N': = Na-N + 2Y^ 7A^e,7 , Ni' = N-2N^ + 2Y^ jN^,^ . (3.23) 

7=1 7=1 

Equation ( p.23| ) shows that the number of holes A"^ in the LWS's II is larger than 
the number of holes in the LWS's I by the additional term 2J2'^=ilNa,'y This 
requires, if we want Eq. ( |3.7| ) for S*" to give the correct eigenvalues, that half of the 
pseudoholes in excess, relatively to the respective LWS I, have /? = +1/2. That is, 
in any LWS/HWS, the numbers of /3 = ±1/2 pseudoholes must be given by 
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N^,p = ^ - P[Na - N] , N^^, = ^ - P[N, - N,] , (3.24) 



the total number of a, pseudoholes being Nj^ = S/j ^a,i3- 

Equation (|3.24|) is a generalization of the pseudohole expressions introduced in 



the context of LWS's I/HWS's I [3- It is convenient to introduce the numbers 
such that 

-| oo 

iV^ = 5" - \Sf\ = - 2/3/«)iV„% - Y: 27iVa,.] • (3.25) 

^ /3 7=1 

This number indicates whether a given eigenstate is an LWS (or HWS). The knowl- 
edge of all the Nj^^^ pseudohole and N^^^ heavy-pseudoparticle (7 > 0) numbers 
provides the knowledge of all pseudoparticle numbers Na^^y (7 = 0, 1, 2, 3...) and 
numbers, the inverse being also true. While the set of Na^^y numbers are directly 
given by the BA solution [Q, the /5-dependent pseudohole numbers Nj^^^ and the 
numbers are not considered by that solution. 

It follows from Eq. (|3.24] ) that the electron numbers are exclusive functions of 
the pseudohole numbers (this had to be so, as required from the form of the 5*0(4) 
generators) and read 

Nr = ^+i: pm, - , n, = ^+y: pk, + . (3.26) 

The pseudomomentum-number operators 

NaAl) = 4,",/3S,«,/3 > = ^l,a,7^g,",7 ) 7 = 1, 2, ... , (3.27) 



and 



iV,,o(g) = l-E^i/3(^)' (7 = 0), (3.28) 

/3 



play a central role in the generalized theory (see Sec. [4.21) . The operator (|3.28|) has 



the following alternative representation in terms of a, pseudoparticle operators 



NaM = blc.A,-fi- (3.29) 

The number operators can be expressed in terms of the pseudomomentum distribu- 
tions as follows 
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Nip = E Km ' ^-,7 = E Km ■ (3.30) 
q <? 

We emphasize that while all Hamiltonian eigenstates are also eigenstates of the 
operator iV^^^ because these numbers are good quantum numbers, this does not 
hold for the operator N^pi^q); obviously, only the states I are eigenstates of the 
latter operator. This fact follows directy from the structure of the state ( |3.18|) . On 
the other hand, all Hamiltonian eigenstates are also eigenstates of the pseudoparticle 
operators Na^.y{q). 

The generalization of the LWS's BA total momentum expression ^ ^, which 
we denote by Pba, to all Hamiltonian eigenstates, requires the addition of an extra 
term associated with rj pairing |Tl], |T^ . This is a direct implication of the commu- 
tators introduced in Sec. (\2-'4) and leads to 



P = PBA + T^lS' + St]. (3.31) 

When this expression provides values of momentum such that \P\ > vr, the total 
momentum is defined as the corresponding value at the first Brillouin zone. In 
operator form and for sub-canonical ensembles such that S*" 7^ 0, and with la given 
by Eq. ( p.l4| ), we have that < N'^ i^, and the present representation leads to 

the following simple expressions for the momentum operator P and its eigenvalue P 



E E QCa,-yNo.M + 7r[iV^"^ + E ^< 

q,a 7=0 ' ^ 7=1 



c,7J ' 



E E QCa,,Na,,iq) + vr E [(1 + 7)iVc,7 + 

q,a 7=0 7=1 



(3.32) 



respectively, where Ca,o = 1 and Cc,7 = —1 for 7 > 0. We note that the momentum 
term, vr J2-y=i [(1 + 7)^c,7 + N^], is always a multiple of ±tt. 

The pseudoparticle perturbation theory introduced in the works [0, 0, |T^, |TB|, 
T7| , |18[ and developed in a suitable operator basis |^, |^, [1^, ^ refers to the 



Hilbert subspace spanned by the Hamiltonian eigenstates I. At finite values of the 
magnetic field and chemical potential and at constant electronic numbers the low- 
energy excitations I are described by pseudoparticle-pseudohole processes relative to 
the canonical-ensemble GS. The latter state, as well as all excited states I with the 
same electron numbers, are simple Slater determinants of pseudoparticle a, levels 
I, I, m, H, HI. That is, in Eq. (CT ) either A^^ 1 = for LWS's or A^^ _i = for 



HWS's. 
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3.4 LWS's and non-LWS's: the criterium 

A LWS is such that 



<i = E7iV.,., (3.33) 
whereas for a S°' HWS we have that 



7=1 



^'-i = E7iVa,7, (3.34) 

' 7=1 

with A^^ 1 > A^'^ 1 for a LWS and A^'^ i > A^^ i for a HWS. 

For the particular case of states I, we have that A'q,^^ = for all 7 = 1,2,... 
branches in the former equations. Therefore, the states I are characterized by zero 
occupancies of the 0,7 heavy-pseudoparticle bands (with 7 > 0). It then follows 
from the general Eqs. ( p.33|) and ( p.34| ) that a S"" LWS I is such that 

<i=0, (3.35) 
whereas for a S°' HWS I, we have that 



iV,^i=0, (3.36) 



Th 

~ 2 

with A^'* 1 > for a LWS I and A^^ 1 > for a HWS I. This leads to the simple 
Slater-determinant form ( p.l8| ) found for the states I . These results also confirm 
that in the case of states I the a,P pseudoholes are such that /? = y but that this 
equality does not hold true in the general case. This is because in the case of the 
states I the Slater determinant ( p.l8|) either involves a,+| or a, — | pseudoholes 



only, whereas in the case of the non-LWS's and non-HWS's or in the case of states 
II (i.e. LWS's or HWS's containing a, 7 heavy pseudoparticles) the state ( |3.18| ) has 
both a, +1 and a, pseudohole occupancies. 

In the case of the non-LWS's and non-HWS's, all or some of the off-diagonal 



generators of the rhs of Eq. ( 3.18|) create new states from a reference BA LWS. The 



3°" non-LWS's and non-HWS's outside the BA are such that 

00 

- Y: 27A^a,7 > K 1 - < 1 1 • (3.37) 

2 2 

7=1 

In contrast, in the case of the LWS's II, the /5 flips generated by these operators 
create a number 7 of (a, ^) pseudoholes for each a, 7 heavy pseudoparticle. Follow- 
ing Eqs. (|3.71 ), (|3.8|) , and (|3.25| ) this is required for the LWS condition 5" = —S^. 
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Therefore, in the latter, cases there is no relation between the quantity la defined 
in Eq. ( |2.14| ) and the pseudohole quantum number /3. 



3.5 The generalized ground state concept 



The possible pseudohole and heavy-pseudoparticle numbers are constrained by Eqs. 
(p.7|), ( |3.8| ), and ( p.25| ) and together with the number of discrete pseudomomentum 
values in each band leads to 4^" possible orthonormal Hamiltonian eigenstates of 
form (|3.18| ), in agreement with the results of Essler et al. |, 0. 

The conservation of the pseudohole and heavy-pseudoparticle numbers permits 
dividing the Hilbert space into subspaces spanned by the set of Hamiltonian eigen- 
states ( 3-181 ) with the same {A*"^^^} and {iV„,^} numbers, which correspond to the 
same sub-canonical ensemble. Obviously, a canonical ensemble (with constant elec- 
tron numbers N^j and thus with constant values of A^'* i — A^'* i for both a = c, s) 
is usually reahzed by several {N^p\, {Na,-^} sub-canonical ensembles. Let us intro- 
duce the generalized ground state (GGS) as the Hamiltonian eigenstate(s) ( |3.18| ) of 
lowest energy in each Hilbert subspace. This type of eigenstate is very useful for the 
construction of the pseudoparticle perturbation theory of Sec. They are of the 
form 



^7(-) „(+) 



\si] n n 



1—1c,,0 l—lpafi 



n n n -'a,cA'. 



q>=ql^> q'=q'^l^ q"=qW,, 



|V),(3.3^ 



and when (i) a = c, 7 > 0, and N^^^ is even or (ii) a = s or a = c and 7 = and 
Na^y is odd (even) and 1°""^ are integers (half-odd integers) the pseudo-Fermi points 
are symmetric and given by 



„(±) 



±[q_ 



C —1 



(3.39) 



where 



(If. 



c,7J 



c,7 



7 > 
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C[Fc,0 — — , 



N, 



a 



<lFs,, = (3.40) 

On the other hand, when (i) a = c, 7 > 0, and N^.^^ is odd or (ii) a = s or a = c 
and 7 = and Na,-y is odd (even) and 1°"'^ are half-odd integers (integers) we have 
that 

Qf2,-/ = iFaa , QfI^^ = -['?Fa,7 " Ca,'y—] , (3-41) 

' a 

or 

?Fa,7 ~ 9F«,7 " Ca,^ — , ?Fa,7 — ^Q'F«,7 ■ (3.42) 

' a 

The GS associated with a given canonical ensemble is always a state I [P], which 
is a particular case of the general GGS expression ( p.38| ). It is useful to denote 
the GS's by \GS] S';, S^) , where and are the eigenvalues of the correspond- 
ing canonical ensemble. The GS expression associated with the [Ic, h) sector of 
Hamiltonian symmetry f/(l) (g) f/(l) reads 

\GS;s:,s:)= n n n n (3-43) 
(-) -{+) ' ' ^ (-) -(+) ' ' ^ 

In all sectors of Hamiltonian symmetry there are states I. In the particular case of 
the SO (4) zero-chemical potential and zero- magnetic field canonical ensemble there 
is only one state I, the pseudohole and heavy-pseudoparticle vacuum, \V), which 
is nothing but the SO (4) GS 0. The same applies to the sectors of Hamiltonian 
symmetry SU{2) (g) U{1) and f/(l) ® SU{2), the GS being always a state I. (In 
addition, in these sectors there is a large number of excited states I.) 

In the case of the two (h) SU{2)(g>U{l) sectors, the GS is both a LWS and HWS 
of the ?7-spin algebra. Therefore, it is empty of c pseudoholes and reads 



g(-) J + ) 



1—1s,0 1—1fs,0 



GS;0,S:)= n n (3-44) 



In the case of the {Q U{1) ® SU{2) sector the GS is both a LWS and a HWS of 
the spin algebra and is empty of s pseudoholes. It reads 
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\GS;S:,0)= li tl i.'JV)- (3-45) 

{-) =(+) ' ' 2 

Finally, the = = (and /i = 0) and = 3^ = SO (4) ground state 
is, at the same time, a LWS and HWS of both the 77-spin and spin algebras p, 
||. Therefore, it is empty of both c and s pseudoholes and it is the vacuum of 
the pseudohole and heavy-pseudoparticle theory. All the remaining eigenstates can 
be described by particular cases of the general expression ( p.l8|) , being the only 
difference the values of the momenta q in the product of operators al^^q^p and 
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Chapter 4 

Pseudoparticle Perturbation 
Theory 



4.1 Pseudoparticle Hamiltonian 

Using the operational representation for the Hilbert space of the one-dimensional 
Hubbard model, we develop a generalized perturbation theory whose small param- 
eter is the density of excited pseudoparticles. A similar theory for the LWS-I has 
been developed previously in the literature 

For simplicity, let us denote the general Hamiltonian eigenstates lip; {Nj^^^}, {Na^^y}) 
by (and the GGS's ( p.38| ) by \GGS)). These states are also eigenstates of the 
a, 7 pseudomomentum-distribution operators (|3.27|) such that 



N^Mli^) = NaMm , (4.1) 

where Na,-y{q) represents the eigenvalue of the operator p.27| ), which is given by 1 
and for occupied and empty values of q, respectively. 



It follows from the form of the state-generator of the rhs of Eq. (|3.18|) that 
all 25*" + 1 Hamiltonian eigenstates constructed from the same S*" = — 5*" LWS by 
applying onto it 1,2,..., 2S°' number of times the off-diagonal operator 5*" (|3.?1 ) have 
the same Na^^^q) pseudomomentum distribution. 

As in the case of the spatial wave functions, the expression of the energy in terms 
of the quantum numbers 1°'''^ of Eq. ( p.ll|) involves the BA rapidities 0, ^ . These are 



complex functions of the numbers If'^- Following our pseudomomentum definition, 
Eq. ( |3.11|) , in our description they are functions of the pseudomomentum q. The 
real part of the BA rapidities, which we denote by Ra^-yiq), are in the present basis 
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eigenvalues of corresponding rapidity operators Ra,'y{q), and obey the eigenvalue 
equation 

R^M\^) = Ra,,{qm , (4.2) 

where for the particular case of Rcfi{q) we also consider the associate rapidity K{q), 
which is defined by the following equation 

R,M = '^^. (4.3) 

u 

The rapidity K{q) is the eigenvalue of the rapidity operator K{q) such that 

K{qM)=KiqM. (4.4) 

The eigenvalues K{q), Rc^^{q), and Rs^^{q) obey the following non-linear algebraic 
equations || 



^^a ^=0 q' J- -t- 7 



7=1 g' 7 

2Re sin-' {u{R^M + ^7)) = ? + ^ E ^c,o(g') tan-^ ^ RcM ' RM) ^ 
+ ^ E E^c,y(?')07,7' ^^(g) - ^c,y(g')) , (4.6) 



and 



2 V-.. , /^._-l/^.,7(?)-^c,o(g') 



E^c,o(g')tan-(^^ 



9' 

1 



+ 7 



where 



7-1 

67,7.(0;) = ey,7(a;)=57,y{2tan-i(|-)+E4tan-i(^)} 
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+ (l-5,„0{2tan-(^^)+2tan-(;^) 



7+7'-|7-7'l ^ 



1 = 1 



+ E 4tan-(^--^)}. (4. 



|7- Yl +2/ 



The limits of the pseudo-Brillouin zones, q^^, associated with the pseudomomentum 
summations are given in Eqs. (|3ll )- (|3T7|) . Ahhough the integral Eqs. ( W^ 



are coupled, note that each of these equations defines the rapidity associated with 
one of the a, 7 bands in terms of other rapidities. 

The key point is that the eigenvalue Na^j{q) is common to the whole tower 
of 2S°' + 1 Hamiltonian eigenstates constructed from the same 5"° LWS. Also the 
associate rapidity eigenvalue Ra^^{q) of Eq. ( ^.2| ) is common to these 25"" + 1 Hamil- 
tonian eigenstates. This follows from the fact that rapidity solutions Ra,-y{q) of Eqs. 
( [4.5|) , (|4.6|) , and are functionals of the Na,^{q) distributions and /3 indepen- 
dent. We thus conclude that the LWS rapidities extracted from the BA provide full 
information on the non-LWS's rapidity operators. 

The rapidity eigenvalues of Eqs. (|4.2|) and (^!^) are independent of the pseudo- 
hole numbers iV^^ and only involve the pseudoparticle distributions Na,^{q). The 
same holds true for the associate rapidity, pseudohole, and pseudoparticle opera- 
tors. Since the SO (4) Hamiltonian (2) and other physical quantities turn out to 
be exclusive functionals of the rapidity operators, in studies involving such quanti- 
ties it is often more convenient to use the 0,7 pseudoparticle representation (with 
7 = 0,1,2,3,...) than the a,P pseudohole and a, 7 heavy-pseudoparticle represen- 
tation (with 7 = 1, 2, 3, ...). 

By direct insertion in each of Eqs. ([4.5|)-(^I7|) of the corresponding q^^^ pseudo- 
momenta we find that the functions i?a^^(g) have for all eigenstates the following 
boundary values at the limits of the pseudo-Brillouin zones 

ii:(g(J)=±7r, i?,,o(gS) = 0, (4.9) 
for a = c and 7 = and 

i?„,^(gW) = ±oo, (4.10) 

for the remaining choices of the quantum numbers 0,7. 

It is useful to consider the rapidity functions K^'^\q) and R'^}^{q) which are the 
GGS eigenvalues such that 
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k{q)\GGS) = K^^\q)\GGS) , (4.11) 

and 

R^,^iq)\GGS) = R'^^l{q)\GGS) . (4.12) 

These are defined by Eqs. ([4.5[ - [4.7|) with the pseudomomentum distribution given 
by its GGS value, which we denote by Nj^}^{q) and is such that 



N^,^{q)\GGS) = Ni%)\GGS) . (4.13) 

Following the GGS expression (iH), N^^^oiq) and N^°^{q) have the simple free 
particle/hole Fermi-like form 



Ni%) = e(Q,(gi+j,,-g)), 0<g<gW 

= e(C«,,(g-glrj,,)), gy<g<0. (4.14) 

For the particular case of a GS, these distributions read 

Nj^^iq) = e(gK,o-9)> 0<g<gl+J 

= ©(g - q'fafi) > giVo < ? < , 
Ni%q) = 0, 7>0. (4.15) 

In Appendix we study the rapidity solutions of Eqs. (^I5|)- (|47| ) both for GGS's 
and for Hamiltonian eigenstates differing from a GGS by a small density of excited 
pseudoparticles. In equations ( [A.?] ) and ( |A.8| ) the parameters Q*-^^ and r*^^ are 
defined combining the two equations 

Q(^) = ir(°)(git!o) , = (gW^) , (4.16) 

with equations (|A.4| ) and ( |A.5|) . For later use, it is convenient to introduce the 



associate parameters 

Q = K(°H?Fc,o) , r^,, = R^SM^n) ■ (4.17) 

In the thermodynamic limit, we can take the pseudomomentum continuum limit 
qj q and the real part of the rapidities become functions of g, which we have 
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called here i?Q,,^(g). In that limit, the BA system of algebraic equations are replaced 
by the system of infinite coupled integral equations presented in Appendix 0. 

We have mentioned that the BA solution is most naturally expressed in the pseu- 
doparticle basis. One consequence of this is the simple expression for the Hamil- 
tonian ( |2.13| ) in that basis. From the SO (4) symmetry we know that the SO (4) 
energy is the same for all eigenstates of a given SO (4) tower. Thus, it suffices to 
know the S0(4) energy for the tower's LWS, which is given by the Bethe ansatz. 
In the pseudoparticle basis, where the non-LWS's/HWS's are characterized by the 
/3 quantum number, the SO (4) Hamiltonian is written, as it should be, in terms of 
the pseudoparticle numbers 



Hso(4) = -2tJ2N,,oiq)cos[Kiq)]+2tJ2N,,^iq)Y.^l-[uiR,M+jt^)]^ 

9 9,7=1 J=±l 

+ ^[^ - - E 7iV,,] . (4.18) 

7=1 

Contrary to the SO (4) term, the chemical potencial and the Zeeman term are f3 
sensitive, because the respective terms in the Hamiltonian ( |2.13 ) do not commute 



with the ladder operators of the SO (4) algebra. Thus, the inclusion in the Hamil- 
toninan ([4.18|) of a chemical potential and a Zeeman terms must be written in terms 
of pseudoholes, and reads 

H = Hso(4) + 2/i E (3N^^p + 2fioH ■ (4.19) 



This is the pseudoparticle expression of the Hamiltonian ( p.l3|) at all energy scales. 



The Hamiltonian ( [4.19| ) then implies that all energy eigenvalues of a given tower 
read E = Eso{i) + 2fi{S'^ + SI) + 2hqH{S'^ + SI). Despite its simple appearance, the 
Hamiltonian ( [4.18| ) describes a many-pseudoparticle problem, the reason being that 
the expression for the rapidity operator in terms of the operators Na^-yiq) contains 
many-pseudoparticle interacting terms. 

Unsurprisingly, it is difficult to solve the BA operator Eqs. ([4.5|), ( |4.6| ), and ( [4.7| ) 



directly and to obtain the explicit expression for the rapidity operators in terms of 
the pseudomomentum distribution operators (|3.27| ). In contrast, it is easier to calcu- 



late their normal-ordered expression in terms of the normal-ordered operators. The 
rapidity operators Ra^-yiq) contain all information about the many-pseudoparticle 
interactions of the quantum-liquid Hamiltonian. There are two fundamental prop- 
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erties which imply tlie central role that the rapidity operators of Eqs. ( |4.11| ) and 
( [4.12| ) have in the present quantum problem: 



(a) As we find in Sec. [4.2| , each of the normal-ordered rapidity operators : -Ra(g) : 
(relative to the suitable GGS or GS) can be written exclusively in terms of the 
pseudomomentum distribution operators ( |3.27|) . 



(b) The normal-ordered version of the SO (4) Hamiltonian (|4.18| ) can be written 



exclusively, in terms of the pseudomomentum distribution operators ( p. 27] ), but all 
the corresponding many-pseudoparticle interaction terms can be written in terms of 
the rapidity operators : Ra,^{q) :. It follows that the rapidity operators commute 
with the Hamiltonian. 

In the ensuing section, we introduce the pseudoparticle perturbation theory 
which leads to the normal-ordered expressions for the Hamiltonian and rapidity 
operators. Despite the non- interacting form of the Hamiltonian eigenstates ( p.l^ ), 
the normal-ordered Hamiltonian includes pseudoparticle interaction terms and is. 



therefore, a many-pseudoparticle operator. However, we find in Sec. [4.2| that these 
pseudoparticle interactions have a pure forward-scattering, zero-momentum trans- 
fer, character. 

4.2 Pseudoparticle perturbation theory 

We will be mostly interested in GS - GGS transitions followed by pseudoparticle 
- pseudohole excitations involving a small density of a, 7 pseudoparticles relative 
to both the GS and GGS distributions. Let us then introduce the normal-ordered 
pseudomomentum distribution operator 

■■N^M :=iV,,,(g)-iVW(g). (4.20) 
When we choose N'^l^^q) to be the GS pseudomomentum distribution (^]T5|) we have 

■.No.M-=NaM^ 7>0. (4.21) 

The normal-ordered distribution ( |4.20| ) obeys the eigenvalue equation 



■.N^M--\^) = ^N^M\^)- (4-22) 
We also consider the normal-ordered rapidity operator 

: Rc^M := Rc^M - R^SM ^ (4-23) 
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where is the GGS eigenvalue of Eq. ( [4.12D 



In the pseudoparticle basis, the normal-ordered rapidity operators : Ra^^{q) : 
contain an infinite number of terms, as we shall demonstrate below. The first of 
these terms is linear in the pseudomomentum distribution operator : Na^^{q) :, 

( [4.20| ), whereas the remaining terms consist of products of two, three, , infinite, 

of these operators. The number of : Na^^{q) '■ operators which appears in these 
products equals the order of the scattering in the corresponding rapidity term. 

A remarkable property is that in the pseudoparticle basis the seemingly "non- 
perturbative" quantum liquids become perturbative: while the two-electron forward 
scattering amplitudes and vertices diverge, the two-pseudoparticle / functions (given 
by Eq. ( |4.5U| ) below) are finite. 

By "perturbative" we also mean here the following: at each point of parameter 
space (canonical ensemble), the excited eigenstates are of form ( |3.18| ) and corre- 



spond to quantum-number configurations involving a density of excited pseudoholes 
and heavy pseudoparticles relative to the GS configuration ( |3.43| ) and ( f4.15| ). We 
introduce the following density 



<. = E<'+ E (4-24) 

a,f} «,7=1 



which is kept small. Here 



and 

1 



ex Ir Eil-<'r(?)]'^^«,/3(g), (4.25) 



n 



a, 7 
ex 



TfT^NaM. (4-26) 



N, 



define the densities of excited a, j3 pseudoholes and a, 7 heavy pseudoparticles, re- 
spectively, associated with the Hamiltonian eigenstate 

In Eq. ( |4.25| ), N^'^p\q) is the GS pseudohole distribution which reads 



<;^(g) 



1 - N^kq) , 
0, 



(4.27) 



if the GS is a 5" LWS and 
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N^'^liq) = 



Klfkq) = l-Ni%), (4.28) 

if the GS is a S" HWS, where Nj^^Q{q) is the a,0 pseudoparticle GS distribu- 
tion (|4.15|) . Our perturbation theory refers to general GS transitions to Hamil- 



tonian eigenstates \ip) involving (a) a GS - GGS transition and (b) a Landau-liquid 
pseudoparticle-pseudohole excitation around the GGS. Therefore, the initial state 
is a GS and we have assumed the GS distribution ( [4.15| ) for the a, 7 heavy pseu- 
doparticles. Following Eq. (|4.15|), we thus have that Na^.y{q) = 6Na^^{q) in the rhs 



ofEq. (1^:2^ ). 

The expectation values of the S0(4) Hamiltonian ( [4.18| ) in the final states {ip) 
are functions of the density of excited pseudoparticles only. This density is given by 



ne. = (4-29) 

a, 7=0 

where n"^^ is given by 



<i' = W^[^- ^^(^)]'^^".o(g) , (4.30) 

for 7 = and by ( [4.26|) otherwise. The density ( [4.29|) is small provided the density 
( [4.24| ) is also small. This then implies that all the densities n'^^ are small and that 
we can expand the expectation values in these densities. The perturbative character 
of the quantum liquid rests on the fact that the evaluation of the expectation values 
up to the n*^ order in the densities ( [4.29 ) requires considering only the corresponding 
operator terms of scattering orders less than or equal to n. This follows from the 
linearity of the density of excited a, 7 pseudoparticles, which are the elementary 
"particles" of the quantum liquid, in 



SN^,,{q) = : N^M : , (4.31) 

and from the form of ([4.26|) and (^4.30|) . The perturbative character of the quantum 
liquid implies, for example, that, to second order in the density of excited pseudopar- 
ticles, the energy involves only one- and two-pseudoparticle Hamiltonian terms, as 
in the case of the quasiparticle terms of a Fermi-liquid energy functional |^ . 
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4.2.1 Pseudoparticle bands and /-functions 

We will develop our perturbation theory considering the two cases where (i) the 
operator normal-ordering is relative to the initial GS and (ii) that ordering refers to 
the GGS associated with a GS - GGS transition. In the following we consider the 
general case of normal ordering relative to a GGS. Later on we will specify when 
that ordering refers either to a GS (which is a particular case of a GGS) or to a 
(non-GS) GGS. 

The expressions for the rapidity operators : Ra,^{q) : corresponds to expanding 
the expressions of the operators : Ra,^{q) : in terms of increasing pseudoparticle 
scattering order. It is convenient to define these expressions in terms of the operators 
: Xa^.^f{q) :. These are related to the rapidity operators as follows 

: Ro^M ■■= K,(l+ ■ ■■) - <7(^) ' (4-32) 

where i?^°^(g) is the suitable Ra^^{q) GGS eigenvalue of Eq. ([4.12|) . 

The operators : Xa^^{q) : contain the same information as the rapidity operators, 
and involve, exclusively, the two-pseudoparticle phase shifts ^a,r,a-',-y'{'ii'i') that we 
define below. In Appendix ^ we introduce ( [4.32|) in the BA equations (|4.5| - |4.7| ) and 
expand in the scattering order. This leads to 

oo 

:X„,,(g):=5:xS(g), 
i=i 

where j gives the scattering order of the operator term X'^)^{q) 
term we find 

^S(g) = |^ E $.,7;a',y(g,g'):iV„',7'(9'):, (4-34) 

where the phase-shift expressions are given below. While the expressions for the 
phase shifts ^a,r,a' ^'{q, q') are specific to each model because they involve the spec- 
tral parameters, the form of the operator term X^l^[q) ( [4.34[ ) is universal and refers 
to all the solvable electronic multicomponent quantum liquids. 
The phase shifts $0,7,0'-^' are given by 

n lis ^ /sin/c sinA;'\ <,r^ 

$c,0;c,0 A;, k') = $,,0;c,0 , , 4.35 

\ u u J 

^c,0;a,y{k,r') = $c,0;a,7 i^^,r'] , (4.36) 



(4.33) 

. For the first-order 
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^a,7;c,o(^ k') = $a,7;c,0 j , (4-37) 

^a,7;a,7' i^, r') = ^a,r,a,Y {r, v') , (4.38) 

where the phase shifts $a,7;«',7' are defined by the integral equations (|A.31|) - (|A.39|) 
of Appendix ^ 

The two-pseudparticle phase shifts can be defined in terms of the phase shifts 

^a,7;a',7' aS foUoWS 

■ (4.39) 

The quantity ^a,-y;a',Y{q, q') represents the shift in the phase of the a', 7' pseudopar- 
ticle of pseudomomentum q' due to a zero-momentum forward-scattering collision 
with the a, 7 pseudoparticle of pseudomomentum q. 

In Appendix 0, we use the Hamiltonian expression ( 4.18| ) in terms of the rapidity 



operators to derive the expression for the normal-ordered Hamiltonian. We find that 
in normal order relative to the suitable GGS (or GS), the Hamiltonian Hso(i)^ Eq. 
( [4.181 ), can be written as 

00 

■.Hso(A) :=E^^'^- (4-40) 

i=i 

For example, the first and second pseudoparticle-scattering order terms read 

^^'^ = E <7(?) ^ ^".7(9) ^ (4-41) 

g,o,7 

and 

^^'^ = ^E E \fc.,r,o.'n'{(l^(l')--KM)----Nc.'A(l')--^ (4.42) 

^^a q,cxa q' ,ai a' 



respectively. All the remaining higher-order operator terms of expressions ( [4. 33] ) 



and (|4.4CI| ) can be obtained by combining the rapidity equations ( ^3|) , (|4.6| ), and 



( [4. 71) and the Hamiltonian expression ([4.18|) . 

In Appendix ^, it is shown that the pseudoparticle bands e|^^(g) can be expressed 
in terms of the phase shifts ( [4.39| ), with the result 



<o(?) = -^-2tcosK^'\q) + 2tJ dk^,,o;c,o {k, K^'\q)) sink , {AA3) 
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-l^iQ) = -lU + 4tRe\/l-u^[R'^>{q)+t'y]^ 

+ 2t _ «,,o;c,7(^,^S(g))sinfc, (4.44) 



and 



e°,(g) = 2t / «,,o;s,,(A:,i?S(g))sinfc. (4.45) 
At the pseudo-Brillouin zones the band expressions have the following values 

<,(9SU) = 0, (4.46) 
and the associate group velocities, Va,^{q), are given by 

. ^ del jq) I. 
= ■ (4-47) 

The two-pseudoparticle forward-scattering phase shifts ^a,r,a',Y{<l, l') defined by Eq. 
( F39| ) and "hght" velocities 

Van = Va,y{qFa,'y) ) (4-48) 

play an important role in the physical quantities when the conformal regime is 
approached 0, ||, 0| • 

All Xj^hq) terms of the rhs of Eq. (|4.33| ) are such that both the / functions 



on the rhs of Eq. (|4.42 ) and all the remaining higher order coefficients associ- 
ated with the operators H^^^ of order j > 1 have universal forms in terms of 
the two-pseudoparticle phase shifts and pseudomomentum derivatives of the bands 
and coefficients of order < j. This follows from the fact that the S*- matrix for j- 
pseudoparticle scattering factorizes into two-pseudoparticle scattering matrices, as 
in the case of the usual BA S'-matrix ||lOl- For example, we show in Appendix 



that although the second-order term H^^'' of Eq. ( [4.42|) involves an integral over the 
second-order function Xj^}^{q) (see Eq. ( |B.12| ) of Appendix ^ ), this function is such 
that H^"^^ can be written exclusively in terms of the first-order functions ( [4.34|) 



+ § E 0{NaJ^ E [X^lijqF.,,)? , (4.49) 
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where \va,'y \ = 6*^,7^0,,^, 9{x) = 1 for a; > 0, and 9{x) = for x < and the "Landau" 
/ functions, fa,r,'y',a'{Q:(l'): are found in that Appendix to have universal form 

00 

+ 27r ^ ^ ^ 0{Na",Y')Ca",Y'Va",Y'^a",Y';a,^{jqFa",Y',q) 
j=±l a" 7"=0 

^a",-y";a',-y'{jqFa",Y',q') , (4.50) 



where the pseudoparticle group velocities are given by Eqs. ( |4.47| ) and ( [4. 481 ). 



We can also write the expression of H^^^ for higher scattering orders j > 2. The 
main feature here is that for all sub-canonical ensembles, energy scales, and pseu- 
doparticle scattering orders, the Landau-liquid terms of that Hamiltonian involve 
only zero-momentum forward-scattering. 

4.2.2 Energy-gap equations 



Let us consider that the normal-ordering of the Hamiltonian ([4.40|) above refers 



to the initial GS. In this case, and in contrast to the low-energy Landau theory 
0, [11], |I2|, |T3|, this Hamiltonian describes finite-energy transitions involving a 
small density of excited pseudoparticles. 

Since the relevant finite-energy excitations involve a topological GS - GGS tran- 
sition (to be discussed in Chapter ^), it is convenient to define the normal ordering 
relative to the suitable final GGS of energy 

uo = Eggs - Egs ■ (4.51) 

In this case the Hamiltonian acts directly on the Hilbert subspace of the final GGS 
where all states have the same choices concerning the integer or half-odd integer 
character of the numbers J°'^ of Eq. ( |3.11| ) . 

Normal-ordering relative to the GGS implies separating the problem (^4.401) into 
(a) a finite-energy term associated with the GS - GGS transitions and (b) a low 
{ijj — ujq) energy theory associated with the Landau-liquid pseudoparticle-pseudohole 
processes around the final GGS. In most cases the excitation energy (^4.51|) associated 
with the GS - GGS transitions is finite. The finite energies are the energy gaps of 
the states II and (or) non-LWS's and non-HWS's relative to the initial GS. 

At finite values of and SI the low-energy Hilbert space is entirely spanned by 



states I [0 |Tl|, |Tj]. We note that, for these states, fixing the electron numbers fixes 
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the pseudohole numbers |TT[ and the a, 7 bands are empty for 7 > 0. Therefore, if we 
fix the electronic numbers we have that at energy scales smaller than the gaps ( [4 .511 ) 
for the states II and non-LWS's and non-HWS's only Landau-liquid excitations (b) 
are allowed. This justifies the Landau-liquid character of the problem at low energies 
|l], ^, |Tl], |16[. Obviously, if we change the electron numbers, there occur at 
low energy GS - GS transitions which are particular cases of the above GS - GGS 
transitions (a). Understanding these GS - GS transitions permits one to express the 
electron operators in terms of the pseudohole operators and topological momentum 
shifts (Chapter H). 

Let us evaluate the general uo expression for GS - GGS transitions such that 
the final sub-canonical ensemble is characterized by vanishing or by small values of 
_i^/Nl^. Here, la = ±1 is defined by Eq. ( |2.14|) . We also assume that both the 



initial GS and final GGS correspond to canonical ensembles such that S*" 7^ and 
belonging the same sector of parameter space, i.e. with the same numbers. We 
note that since GS's are always S"" LWS's I or HWS's I there are neither a, — y 
pseudoholes [see Eqs. ( p.35| ) and ( |3.36| )] nor a, 7 heavy pseudoparticles (7 > 0) 
in the initial GS, and that = [see Eq. (|3.25|) ] for that state. From Eqs. 
(CT)-(ra), we find for the GS - GGS gap 

^0 = E e%(gFa,o)AiV,,o + E <7(0)iVs,7 + 2fiASt + 2fioHAS: . (4.52) 

a 7=1 



Inserting in Eqs. ( |3.24]) both the pseudohole expressions (|3.23| ) and the expressions 



S^^ = -l[Na~ N] and = -|[iV| - N^], we find, after use of Eq. (pl), the 
following general expressions for AA^^o and AA^^^o 

7=1 

- - AiV^\ ] - 2iV,^ - E 27iVc,7 , (4.53) 



2 ' 2 , 

7=1 



and 



AN. 



s,0 



Y^ILAS: + K]-Y: 7iVc,7 - E(l + 7)iV.,7 

a 7=1 7=1 

^E[^i-i--^^'i- + 2A^«] 

2 ^ °' 2 "'2 

E7iVc,7-E(l + 7)iVs,7- 

7=1 7=1 



(4.54) 
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Using these and Eq. ( p.25D in the rhs of Eq. ( [4 .521 ) leads to the gap expression 
(pjm of Appendix M. 



On the other hand, generalization of the results contained in the work ||T2[ for the 
chemical-potential and magnetic-field expressions to all sectors of parameter space 
leads to 



2fxo\H\ 



Fs,0 



(4.55) 



From the use of Eqs. ( [4.53|) and ([4.54|) in the gap equation ( |4.52|) , we can rewrite 
Eq. (^l0| ) as follows 



uo = 2|/i|iV^';_^ + 2fio\H\[N^^_^ + iV,,,] + <7(0)iV,,^ , (4.56) 

7=1 7=1 

which is the general expression for the gap for GS - GGS transitions. It follows from 
Eqs. ( p. 251 ) and ( |3.7| ) that the gap expression (|4.56| ) can be rewritten as 



^0 = 2|/i| 7iVc,7 + N:)+2f,o\H\ (E(l + 7)iVs,7 + E <,ms,, • (4.57) 

7=1 7=1 7=1 



For GS - GGS transitions to pure LWS's or HWS's II, the use of Eq. leads 



to 



K = 0, (4.58) 

for both a = c,s. It follows that for such transitions the gap is given by expres- 
sion ( [4. 571) with = and A^^ = 0. Therefore, the single (7 > 0) 0,7 heavy- 
pseudoparticle gap is 2'y\fi\ and 2(1 + 'y)fio\H\ + e°^(0) for c and s, respectively. 
Given a final sub-canonical ensemble we can write 



: H := Ho + Hiandau , (4.59) 

where Hq has eigenvalue ojq and corresponds to the GS - GGS transition (a) and 
H Landau is uormal-ordcrcd relative to the GGS and is of the form 

HLandau = + , (4.60) 

where 

H^^= J2 ■■ ■■ ^ (4-61) 

g,o,7 
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e.M = ^aliq) - e^aii<lFa,o) , = egUg) - 5„,.6gUo) , (4.62) 



and H^'^^ is given in Eq. ( 4.42| ). The Hamiltonian ( [4.60| ) describes the Landau-hquid 



excitations (b) relative to the GGS. Here ( [4.60| ) are the Landau-hquid Hamiltonian 
terms which are relevant at low energy {uj — ujq). Therefore, the corresponding 
second-order Hamiltonian is suitable to study the physics at low positive energy 
above the gap, i.e., small {u — luq). 

In general, the different final GGS's of the Hilbert subspace where a given initial 
GS is transformed upon excitations involving a small density of pseudoparticles have 
different energies ( |4.56| )-( p:.57| ). This implies that the study of the quantum-liquid 
physics at energy-scale ujq involves, in general, transitions to one sub-cannonical 
ensemble only. However, if two or several possible final GGS's had the same energy 
gap ujq, the physics would involve the Hilbert subspace spanned by all Hamiltonian 
eigenstates associated with the corresponding different sub-canonical ensembles. 

For sub-canonical ensembles such that S*" 7^ and characterized by small values 
of __ia/Na (where is given in Eq. (BTTl)), the gap expression ( [4.56| ) leads 
directly to 

00 

^0 = 2|/i|iV^^_^ + 2/io|/f| [iV^"_^ + E ^^.7] + E 4°^(0)iV,,^ . (4.63) 

7=1 7=1 



There is a remarkable similarity between the general Hamiltonian ( [4. 591) and 



its version in the Hilbert subspace spanned by the states I [11], Moreover 



in the case that the GGS is the GS it self, the gap (|4.56|) -( l4757| ) vanishes and the 
Hamiltonian (|4.59|) reduces to the above low-energy Hamiltonian. In this case it 
refers to Landau-liquid a,0 pseudoparticle-pseudohole excitations around the GS. 
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Appendix A 

Normal-ordered rapidity operator 
expressions 



Following the discussion of Sec. \l.2[ the perturbative character of the system implies 
the equivalence between expanding in the pseudoparticle scattering order and/or in 
the pseudomomentum deviations ( [4.31 ). In this Appendix we give a short descrip- 



tion of the calculation of the normal-ordered operator expansion for the pseudopar- 
ticle rapidities ( [4.33| ). We focus our study on the evaluation of the functions X^^^(g) 



and Xl^}^(q) and associate functions Qajyil) and Q^^^(g) (eigenvalues of the operators 
X2-\q) and XP{q) and Q^^\q) and Q^^Kq), respectively). 

We evaluate here the first-order and second-order terms of the eigenvalues of 
the operator (|4.33| ). Equation (|4.22|) then allows the straightforward calculation of 



the corresponding operator expressions. All the quantities are given in the limit 
established by Eq. 

In the thermodynamic limit, Eqs. (|4.5|)-(|4.7|) lead to the following equations 



K{q) = q-- / dq'Ns^oiq') t^n'^R^M ' Rs,o{q')) (A.l) 



q = 2i?esin ^ (u{Rc,-y{q) + i"y 
1 



dqN^fl{q)ia.n '- 



(A.2) 



and 
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Normal- ordered rapidity 



TT 



.(+) 



9c, 



1 + 7 



+ 



1 rqi'\^ III I \ 



(A.3) 



In what follows, it is useful to define the GGS rapidities K^^\q) and R'^l^{q) in 
terms of their inverse functions 2'7rpQ, -y(r) as 



dk2npcfi{k) , 



(A.4) 



(A.5) 



We start by considering the GGS eigenstate rapidities of Eqs. ( [4.11| ) and ( [4.12| ). 
If we insert in Eqs. ([4.5p-(^77[) the GGS distributions ( [4.14D , after some algebra we 
find that the functions of the rhs of Eqs. ( |A.4| ) and ( [A.5| ) defining the inverse of 
these GGS rapidities are solutions of the following integral equations 



2vrpe,o(A;) = 1 + 



COS K f^s,o 



U 



' dr- 



\ -\. (^f — sinfcp 



(A.6) 



27rp,,^(r) = 2i?e( 



u 



1 — u^[r + 



Q(-) 



dk- 



2Txpcfi{k) 



TT'-f 



1 + 



(A.7) 



and 



27ips 



where 



dk- 



2-Kpcfi{k) 



vr(l + 7) 



1+7 



f ^eSV7(---')2vrp,o(r') 

(A.l 



dx 



7-1 
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S (i7-Yi+20[i+(p^;f=H^' ^^-'^ 

is the derivative of the function ( [4.8| ) and the parameters Q*-^-* and r^^^ are defined 
combining Eqs. (^), (pTSl) , and (|4.16|) . 

Let us now consider small deviations from a GGS or GS. The eigenvalue form of 
Eq. ( ^ ) is 

SK{q)=K^'\q + 6XM)-K^'\Q)^ (A-IO) 
for a = c and 7 = and 

5R^,,{q) = <^(g + - , (A.ll) 

for all remaining values of the quantum numbers a and 7. Here 6K{q), SRa^-yiq), 
and are the eigenvalues of the operators : K{q) :, : Ra,'y{q) '■, and : Xa^-yil) '■, 

respectively. From Eq. ( |4.33| ) 5Xq, ,y(g) can be written as 

= X« (g) + X^^) (g) + ... , (A.12) 

where Xj^^^{q) is the eigenvalue of the operator Xj^^^{q). Expanding the 6Ra,^{q) 
expressions ( [A.lOp and ( |A.11|) we find 



K{q) = J2K^Hq), (A.13) 

i=0 

and 

00 

i?„,,(g) = Ei?«(g), (A.14) 

i=0 

respectively, [and SK{q) = EZiK^'Kq) and 5Ra,^{q) = where the 

zero-order GGS functions K^'^\q) and Ra^-^iq) are defined by Eqs. Eqs. ( |A.4| ) 



and ( [A.5|) . From the resulting equations we can obtain all derivatives of the GGS 
functions K^'^\q) and with respect to q. The terms of the rhs of Eqs. (|A.13D 



and ( |A.14| ) involve these derivatives. For instance, the first-order and second-order 
terms read 

A-'"(,) = ^<;Jto). (A.15) 
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(g) = ^^xS(g) , (A.16) 



and 



flW (,) = + l^^lM[A-a)(,)]^ , (A.18) 

respectively, and involve the first and second derivatives. 

From Eqs. (^)-(^) [with N^^^{q') given by the GGS distribution (75)] we 
find that the first derivatives '^^'^ ^'^^ and '^•^"■t'-'^^ can be expressed in terms of the 

aq aq ^ 

functions as follows 

and 



(A.19) 



dR^^Uq) 



dq 27rp,,^(i?i°^(g)) 



(A.20) 



The second derivatives ^^g^^'^'^ and ^ ^"'i ^'^^ then read 



dq^ ~ "[27rp,,o(i^(°)(9))]^ rf^''^""'"'^"^ ' ^ ^ 

and 



dq' [27rp,,,(/?a(g))]3 dr 

By introducing both the distributions 



N^M = Nj^hq) + 5N^,,iq), (A.23) 



and the first-order and second-order functions ( [A.15D -( [A.22D into Eqs. ( |A.1D -( |A7^ 



we find after expanding to second order that for j = 1 and j = 2 the functions 
X^}^{q) can be written as 
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where Y^^^^{q) is even in q 

and does not contribute to the physical quantities to second scattering order (and 
to second order in the density of excited pseudoparticles) . Up to this order only 
Qa^'yil) contributes. 

To derive this result we have expanded the expression for the even functions 
^i/] (g) for j = I and j = 2 to second order in the density of excited pseudoparticles. 
Following the perturbative character of the quantum liquid in the pseudoparticle 
basis, the obtained expression is exact up to j = 2 pseudoparticle scattering order. 
Moreover, we find that the even function Y^^^ (q) does not contribute to the physical 
quantities up to that order and can, therefore, be omitted. 

Let us introduce the functions Qajyii^) such that 

QS(^) = '5S(^S(?))' (A-26) 
for all values of a and 7, with i?^°Q(g) = EI}J^Ll^sl^ it is also useful to define the 
function Q^^\k) such that 

Q''\k) = Q^IC^) . (A.27) 
u 

By introducing both the distributions ( |A.23|) and the first-order functions defined 
in Eq. ( [A.26| ) into Eqs. ( [A.1| )-( [AI3| ), we find after expanding to first order that the 
functions Qa)yii^) are defined by the following system of coupled integral equations 




7'=i '^c,y 
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and 



,.(+) 

{-) 
c,0 



dr'- 



TT'-f 



1 + 



' r—r' \ 2 



(A.29) 



(i)( 



dq5Ncfi{q)— tan" 



+ E 



7'=0 



•5,7 

(-) 



+ 



,(+) 



(-) 

c,0 

{-) 
s,0 



/TT 



1 + 7 



vr(l + 7) l + 



/TT 



(A.30) 



(The use of Eq. ( [A.28| ) in Eqs. ( [A.29| ) and ( |A.3CI| ) allows the expression of both 



gW(^) and Q«( r) in terms of free terms and integrals involving Q 



(1), 



Com- 



bining Eqs. ([OsD - dOoD with Eqs. ( [OeD and ([A^TD leads to Eq. (|]3|) with 
the phase shifts defined below. Note that at first order we can either consider the 
function X^'Jyi.Q) '^^ the associate function Q'^ajyil) of Eq. ( |A.24|) . Both functions are 
of the form ( [4.34|) . However, the general expression for the phase shifts associated 
with the function Xj^l^{q) of Eq. ( |4.34| ) have extra terms. These arise from the func- 
tion Y^^^{q) of the rhs of Eq. (|A.24|) . If we expand the physical quantities involving 
the phase shifts to first (and second) order in the density of excited pseudoparticles 
these extra terms lead to vanishing contributions. For instance, expression ( |4.34| ) is 
identical if we use in it either choice for the phase shift expression. For simplicity, 
we omit here the phase-shift extra terms associated with the function Y^^^{q). We 
find that the phase shifts associated with the function Q^ajyil) defined by the 
following coupled integral equations 



*c,0;c,o (r, r') = - j 

TT J r 



S+) 



J-) ]_ _|_ _ ^"y 



(A.31) 



tan 

TT 



-1 



7 



.(+) 



- „^sfi;c,^{r",r' 



{-) 

s,0 



dr" 



1 + (r — r 



f/\2 



(A.32) 



c,0;s,7 



(r, r') 



tan 

TT 



1+7' 



.(+) 



. „^sfl;s,^{r",r') 



{-) 



dr" 



1 + r 



(A.33) 
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^c,r,cAry) = ^tan-i(^) - ± j dr" , (A.34) 



1 




TT , 




1 




TT 





1 ri-!-J , „$^,.v(r",r'^ 
y[l- 

,..(+) 



*s,7;c,o l?^, ) = tan (— — ) + -/,, ar 



rfr"|..,o;c,o (r", r') ^^^^^^^ Z , (A.37) 

s,0 



(l + 7)[l + (^)^] 



,(+) 

3,0 



;,0 



,{+) 

s,0 



®l+7,l+7'(^-^') , 1 ri^.' „ <le,0;s;7'(^",O 



1 


rc,o 


TT 









For 7 = Eqs. ( |A.371 ) and (|A.39|) can be rewritten as 



(A.39) 



1 1 r^^^ 

^'.,0;c,o(r,r') = — tan-i(r-r') + - /)^.,0;c,o x') , (A.40) 



and 



- 1 _i,r-r' 1 riV , ,Mn-Hr" - r' 

^s,0;s,oir,r') = -tan (^^)-^/,_, 



vr ' 2 ' 7r2 7r(-) 1 + (r - r")^ 

+ tir"G(r,r")<l.,o;s,o(r",r') , (A.41) 
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and the kernel G{r, r') reads 



2tt 



1 



1 + ((r-r')/2)2 



1- o {t{r)+t{r')+ ^ ' ^ ' 



, (A.42) 



with 



and 



t{r) = - 

IT 



tan"^(r + ri'^}) - tan"^(r + rj^o^' 



vr 



/(r) = - ln(l + {r + r[:\)^Y) - ln(l + (r + r^o 



(A.43) 



(A.44) 



In order to evaluate the second-order functions Qajy{q) of the rhs of Eq. (|A.24|) 
for j = 2 we introduce the functions ( [A.15| )-( [A.18| ) and distributions ( |A.23| ) in Eqs. 
( |A.l| )-( |Al^ ). Expanding to second order we find after some algebra 



1 d 



where 



It is also useful to define the function Q'^'^'*\k) such that 

, sin k , 



u 



(A.45) 



(A.46) 



(A.47) 



The functions Q^^^"^ (r) are defined by the following system of coupled integral equa- 
tions 



27rp,,^(i?i°](g)) 7r(l + 7) 

/(_, dq6NcM 



7=0 5s, V 
,(+) 



1 + (IZ^)2 



7=1 " yc,7 



+ E 



E 



27rpc,7(^c°7(g); 7r7 

[Qi!^(jgFs,7)] 



1 + 



+ E 



27rp,,7(r,,7) 27r(l + 7) [l + 
1 



'1 27rpc,7('^c,7) i=±i 2vr7 



1 + 
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+ 



.(+) 

"° dr' 

TT [1 + (r - r 



A21 



(A.48) 



and 



9c,o , , .COS 



7r7 



1 + 



7' = 1 

cosQ 1 

u 



E 



27rpc,o(Q) j=±i 27r7 



7'=1 
.(+) 



Q 



7r7 



(A.49) 



s,7 



,(+) 
'c,0 



dq6N,M 



cosi^W(g) Q2(g) 



« 27rpe,o(/^(°H9)j7r(l + 7) 



1 + 



. 1+7 ) 



+ 



E 

7'=0-"? 

cosQ 



•5,7' 

(-) 



u 



2vrpc,o(Q) iSi 
1 1 



27r(l + 7) 



+ 



^0 27rp,,y(rs,7') 2 ^-^ 



1 



--±1 



c,0 
(+) 



I 1+7 J 

/;:f^r'Q&^(r')^eW,,(r-r'). 



(A.50) 



As for the first-order case, the use of either the function Q^^^(g) defined by Eqs. 
(|05|)-( |X:i^ ) or of the full function Xj^)^{q) of Eq. ( [OlD for j = 2 leads to the 
same results for the physical quantities up to second order in the density of excited 
pseudoparticles. 
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Note that the free terms of Eqs. (|A.48| )-( [A.50| ) involve the first-order functions 
only. This implies that the unique solutions of these integral equations can be ex- 



pressed in terms of the first-order functions. Therefore, following Eqs. (|4.34| ) and 
( A. 24 ) the second-order functions can also be expressed in terms of the pseudopar- 
ticle phase shifts. 



Appendix B 

Normal-ordered Hamiltonian 
expression 



In order to derive the first-order and second-order Hamiltonian terms of Eqs. ( [4.41 
and ( [4.42| ), we again consider eigenvalues and deviations. The energy associated 
with the Hamiltonian ( [4.18| ) reads 



Eso{A) = — ^ y,.) N^,o{q) cos K{q) 

1 ^, 



27r ^^i-'gc.v 

+ f/[x-^-E7iVc,.]. (B.l^ 

^ ^ 7=1 



To caculate the bands ( |4.43| )- (|4.45| ) we introduce in this energy Eqs. ( |A.15|) , ( |A.16|) , 



and ([A.23|) and expand the obtained expression to first order in the deviations with 



the result 



a4o(4) = ^{/f^^c,o(g)[-2tcosKW(g)-|] 

7=1 •'Icri 

2t / dkQ^^\k)smk 

^( 



/q(-) 

„{+) 



7=1 -^^-r \/l-n2r + «7F 
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where the functions Q^^\k) and Q^^}^{r) are defined by Eqs. ( [A.27| )-( |A.29| ). The use 
of Eqs. (lOl , (pOGi) , ( [/OTD , and in (^) leads after some straightforward 

algebra to 



5: eO^^(g)5iV,,,(g) 



(B.3) 



with the bands given by Eqs. (|4.43| )- (|4.45| ) in terms of the phase shifts ( |A.31|) -( |A.39| ). 
An equivalent representation for the band energies is 



e%iq) = -| - 2tcosir(°)(g) - J^^ rfr2tr/,,o(r)^ tan'^ (r - i?S(g)) , (B.4) 



and 



-7f/ + 4tRe\/l - u^[R§i{q) + i-f]^ 



/ dk— tan 

/q(-) 71 



1 



7 



(B.5) 



Q(+) -j^ /sinfc 

-2t / rfA;-tan~^ 

Jq(-) tt \ 1 + 7 

.(+) 



sin /c 



.(+) tan 
" dr'- 



-1 (r^^<M 

1+7 



where the functions 2trjcfi{k) and 2tria^^{r) are defined by the integral equations 

2t?7s,o('^) 



(B.6) 



/,x „ . , cosfc /■ 
2tr]co{k) = 2tsmfc H / 



COS f ^s,o 

, dr- 



1 _|_ _ sinfc\2 



B.7) 



-AtRe(^ 



V? [r + ^7] 



+ 



:>(+) 



2tr]cfi{k) 



:b.8) 



Q(-) 



7r7 



1 + 



and 
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2t / dk- 



sin k 



vr(l + 7) 



1 + 



1+7 



,{+) 



;:;dr'2tr/.,o(r')^[eg+,(r'-r) 



s,0 
J+) 

(-) 



Co' 7r2(l+7) l + [l + (r'-r")2] 



(B.9) 



The gap expressions ( |4.56| ) and ( ^4.57|) can be rewritten in terms of pseudoparticle 
bands as 



and 



7=1 



(B.IO) 



[KMfc,,) + e° o(gFs,o)] (E 7iVc,7 + K) 

7=1 

<o(9Fs,o)(E[l + 7]iV.,7 + E <7(0)^' 



s,7 1 



(B.ll) 

7=1 7=1 

respectively. In order to derive the expression for the second-order Hamiltonian 
( [4.42| ) and associate / functions ( [4. 501) we expand the energy ( p.l|) to second-order 
with the result 



^-^50(4) 



+ 



+ 



§{/:5<'«*'v,„(,) 



27rp,,o(i^(°Hg) 



E /(-)' dq5N,,^{q) 
2tsing 1 ^ (1) , 2 

^^^2^E[^^^°^^^-^°)] 



l-M2[i?(°}(g)+^7]^ 



^^[rc,^ + ^7] 



Q(-) 



dkQ^^^*\k)2t sin A; 
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Normal-ordered Hamiltonian 



E 

7=1 



,(+) 



(-) 

c,'y 



[r + Z7] 



}• 



(B.12) 



Inserting the suitable functions in the rhs of Eq. ( B.12|) , performing some integra- 
tions by using symmetry properties of the kernels of the integral equation (|A.28|) - 
( |A.3CI| ) and ( [A.48| )-( [A.50| ), and replacing deviations by pseudomomentum normal- 
ordered operators ( [4.20|) we find after some algebra expression (|4.49|) . Note that 
replacing in Eq. ( [4. 491) the function X2')y{q) by the associate function Q^}y{q) leads 
to the same result. By the use of Eq. ( |4.34| ) expression ( |4.42| ) can be rewritten in 
terms of the / functions ([4.50|) as given in the rhs of Eq. ([4.42|) . 
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Chapter 5 



Excitations and Finite-Size Effects 



5.1 Introduction 

In this Chapter we introduce the quasiparticle concept for the ID Hubbard Model. 
Exception made to Sees. |^ and p75| , we deal in this Chapter with LWS's I/HWS's 



I only. As we shall see, these quasiparticles have not exactly the same physical 
meaning as in Fermi-liquid theory |l|, ^, || . In the latter theory, the quasiparticle 
concept is used to describe low-energy excitations at constant or variable electronic 
and magnetization densities, n and m, respectively. The quasiparticle in a Fermi- 
liquid is related to the electronic operators trought a factor Z -the renormalization 
factor. The value of Z is determined by the electronic correlations. 

In our Landau-liquid description of the Hubbard chain the quasiparticles do not 
describe the elementary excitations at constant n and m. These are described by 
the a pseudoholes and by the a, 7 heavy-pseudoparticles, introduced in Chapter H. 
In contrast, excitations that change n and m are described by composite objects 
of pseudoholes, that we call (in this restricted Fermi-liquid sense) quasiparticles. 
Due to the non-perturbative character of the one-dimensional many-body problem 
IP, 1^, H, the renormalization factor that connects the quasiparticle and electronic 
operators vanishes as the Fermi points are approached 0. 

The quasiparticle operator is defined as the generator of a ground-state-ground- 
state transition involving the addition or removal of one electron to or from the 
many-body ground state. We find that the quasiparticle momentum equals the 
expected electronic momentum k = ±kp„. It is the combined momenta (i) from the 
addition or removal of pseudoparticles and (ii) from the global pseudomomentum 
shifts of ±Tr/Na which leads to this final result. This type of pseudomomentum shifts 
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are also present in excitations connecting the ground state and LWS's II (or HWS's 
II). We will illustrate the importance of these momentum shifts ||T0| , pJ| in the limits 
n = 1 and U/At^ ^ 1, where the model is closely related to the anti-ferromagnetic 
Heisenberg chain. 



5.2 Ground-state — ground-state topological exci- 
tations 

The quasiparticle operator, c^^^^^, creates one quasiparticle with spin projection a 
and momentum kpa- As in Fermi liquid theory, such operator is the generator of 
the following ground-state-ground-state transition 

cL.JGS; iV., iV_.) = \GS- + 1, N_,) . (5.1) 

The quasiparticle operator defines a one-to-one correspondence between the ad- 
dition of one electron to the system and the creation of one quasiparticle. The 
momentum of the ground states \GS\ N^-, N^a) associated with canonical ensembles 
of U(1)®U(1) symmetryis given in Table p?T| . Since we are studying ground-state- 
ground-state transitions, a, 7 > 1 heavy-pseudoparticles are not involved. 





(-1,-1) 


(-1,1) 


(1,-1) 


(1,1) 


(A)P 


±2kF 


±2kF 


±[2n - 2k f] 


±[2n - 2kF] 


(B)P 


±.kp^ 


±kFi 


±kFi 


±kFi 


(C)P 


±kFi 


±kFi 


±kFi 


±kFi 


(D)P 















Table 5.1: Values of the ground-state momentum in the four (/, /') sectors of Hamil- 
tonian symmetry U {1)®U{1). The different momentum values correspond to the fol- 
lowing parities of the numbers Ns^ and Ng^, respectively: (A) even, even; (B) even, 
odd; (C) odd, even; and (D) odd, odd. In the case of the (1, ±1) sectors, if kpa > vr 
then ±kFa should be replaced by the first-Brillouin-zone momenta ±[2tt — kpa]- 

The study of the ground-state momentum ( |3.32| ) confirms that the relative mo- 
mentum of ground states differing in the number of a electrons by one equals the 
U = Fermi points, ie AP = ±kFa- 
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It is interesting to study the expression of the a quasiparticle and quasihole 
operators in the a, (3 pseudohole basis for all sectors of Hamiltonian symmetry. 
Since we are discussing the problem of addition or removal of one particle, the 
boundary conditions ( p.l2| ) and ( |3.13| ) play a crucial role. When we add or remove 
one electron from the many-body system we have to consider the transitions between 
states with integer and half- integer quantum numbers /°'°. The transition between 
two ground states differing in the number of electrons by one is then associated with 
two different processes: a backflow in the Hilbert space of the a, (3 pseudoholes with 
a shift of all the pseudo momenta by and the creation and (or) annihilation of 
one pair of c and s pseudoholes at the pseudo-Fermi points (or at the limit of the 
pseudo-Brillouin zone for the s pseudohole). 

The backflow associated with a shift of all the pseudomomenta by is de- 
scribed by a topological unitary operator such that 

Obviously, the pseudohole vacuum \V) is invariant under this operator, this is 
V^|l^) = \V). It is simple to find that reads 

= exp |-^4±^^^^^ag,„,^| . (5.3) 

The structure of this operator reveals that the pseudomomentum of all pseudoholes 
are translated by ±11/ Na- Adding all pseudohole contributions gives a large mo- 
mentum. This large-momentum excitation induced by the operator (|5.3| ) is the a 
topological momentum shift. Note that in the present case of the Hilbert subspace 
spanned by the states I of the (/, /') sector only the value /3 = | for a = c and the 
value /3 = I for a = s contributes to the (3 summation of Eq. (|5.3|) 

In addition to the topological momentum shift, the quasiparticle or quasihole 
excitation includes creation and (or) annihilation of pseudoholes. The changes in 
the pseudohole and pseudoparticle numbers and the corresponding changes in the 
values of S'^, S% S'^, and are given in Tables |5]^ and ^]3| for the ground-state 

- ground-state transitions (iVt^^i) ^ (^T ± ^'^i) (^T'^i) i^h^l ± 1)' 
respectively. 

We consider below the expressions for the quasiparticles c|.^^ ^ and c|,^^ associ- 
ated with the transitions (A^t^^i) ^ (^T + l^^i) and (A^t^^i) (^T^^i + 1)^ 
respectively, and the quasiholes Ckp^^i and Ckp^^i associated with the transitions 
(A^t, Ni) ^ {N^ -l,Ni) and {N^, N^) ^ {N^, - 1), respectively. 
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(-1,-1) 


(-1,1) 


(1,-1) 


(1,1) 




Tl 


Tl 


±1 


±1 


AN, 


±1 


±1 


±1 


±1 


AN* 














AN^ 


±1 


±1 


±1 


±1 


ANs 





±1 


±1 





an; 


±1 








±1 


AS" 


Tl/2 


±1/2 


±1/2 


±1/2 


AS"^ 


±1/2 


±1/2 


±1/2 


±1/2 


AS' 


±1/2 


±1/2 


±1/2 


±1/2 


A5| 


Tl/2 


±1/2 


±1/2 


±1/2 



Table 5.2: Changes in the numbers of pseudoholes, pseudoparticlcs, pscudoparticle 
orbitals, and in the values of S'^, S^, S', and S^ in the ground-state-ground-state 
transition (7V|, N^) {N^±l, 

We emphasize that because the initial ground state for the above two quasipar- 
ticles and two quasiholes is the same, the a quasiparticle and a quasihole momenta 
differ by ±;^. Therefore, the corresponding quasiparticle and quasihole expressions 
are not related by an adjunt transformation. On the other hand, the operators 
^±kFa,<T ^±kFa,T assoclatcd with the transitions {N^, N^^) — > {N^ ± 1, N^^) and 
(A^o-±l, N-a) — {^a: ^-a) ^rc obviously related by such transformation. In this case 
the initial (final) ground state of the electrons (holes) is the final (initial) ground 
state of the holes (electrons). Moreover, let us consider the set of four operators 
^±kF^,V ^ifeFT.T, ^±feFi,i' ^i^Fi-i such that the creation operators act on the same 
initial ground state (iV-f, iV|) transforming it in the ground states {N-^ ± 1, iVj and 
{N^,N^ ± 1), respectively, and the hole operators act on the corresponding latter 
states giving rise to the original ground state. Let us consider the reduced Hilbert 
subspace spanned by these three ground states. If we combine that with the elec- 
tron and hole expressions introduced below, it is easy to show that the corresponding 
quasiparticle and quasihole operators cl_i^^^^p c±kp-f,i, and c±kpi,i obey the 

usual anticommutation relations. 

The two electrons and two holes refer to the same initial ground state. The 
pseudo-Fermi points and pseudohole- Fermi points of the expressions below refer to 
that initial ground state. On the other hand, s pseudohole creation and annihilation 
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(-1,-1) 


(-1,1) 


(1,-1) 


(1,1) 






Tl 


±1 


±1 




±1 


±1 


±1 


±1 


AN* 














AN^ 


Tl 


±1 


±1 


±1 


ANs 


±1 








±1 


an; 





±1 


±1 





AS'' 


Tl/2 


±1/2 


±1/2 


±1/2 


ASt 


±1/2 


±1/2 


±1/2 


±1/2 


AS' 


Tl/2 


±1/2 


±1/2 


±1/2 


ASI 


±1/2 


±1/2 


±1/2 


±1/2 



Table 5.3: Changes in the numbers of pseudoholes, pseudoparticles, pseudoparticlc 
orbitals, and in the values of S'^, S^, S', and S^ in the ground-state - ground-state 
transition (A^^, N^) (iVf , A^; ± 1). 

operators at the limits of the pscudo-Brillouin zones refer to the final and initial 
ground states, respectively. In the case of the (/, /') sectors of Hamiltonian symmetry 
U{1) U{1) we consider that the initial and final ground states belong the same 
sector of parameter space. In the case of the (/') sectors of Hamiltonian symmetry 
SU{2)®U{1) [or (/) sectors of Hamiltonian symmetry U{l)®SU{2)] we consider that 
the initial and final ground states belong to sectors of parameter space characterized 
by the same value of (/') [or (/)]. We present below the electron and hole expressions 
found for different initial ground states in the nine sectors of parameter space. It is 
useful to introduce the pseudohole- Fermi point q^^^ such that = ±27i /Na 
For initial ground states in the (—1,-1) sector of Hamiltonian symmetry U{1)^ 
U{1) we find 

^I'l 'iPc 2 Qs 2 IFc 2 ^Fs 2 

for the electrons and 

^Fc '^'2 IS , , 2 Vpc ' ' 2 ^Fs ' ' 2 



for the holes. 
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For the (—1, 1) sector we find 



c\u^^^ = V^a.{±) ia_(±) 1, c+L,^, I = a_(±) xV^a\^. (5.6) 
for tlie electrons and 



Cifc^T.T = K^«l(±) , _ 1 « J±) , 1 > Cifc^ia = « J±) , _ 1 {±) , 1 > (5-7) 

2 ^Fs '*'2 ^Fc ''-'2 '^s ' '2 

for the holes. 

In the (1,-1) sector the result is 



'iPc '^'2 ^Fs 2 ^Fc '^'2 1 > 2 

for the electrons and 



^Fc ' '2 ^Fs ' ' 2 ^Fc ' '2 i^; 2 

for the holes. 

The expressions for the (1, 1) sector are 



for the electrons and 



c±kp^,1-a'g{±)^iV^^aj±) 1, Cifcj, I - (±) lO (±) 1 , (5.11) 

' IFc '^'2 ' '2 ^Fc ''-'2 ^Fs 2 

for the holes. 

According to Eqs. (|5^ )- (|5.11| ) the a quasiparticles and quasiholes are many- 
pseudohole objects which recombine the colors c and s (charge and spin in the limit 
m = n-^ — 111 —>■ 0) giving rise to spin projection | and | and have Fermi surfaces at 

±kFa- 

Similar expressions can be derived for the sectors of parameter space where the 



Hamiltonian (|4.19| ) has higher symmetry. We start by considering the sectors of 



Hamiltonian symmetry SU{2) (E) U{1) where 
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= -q'fj = = -Qi-^ = ^[1 - ■ (5.12) 



•a 



For ground states of the I' — —1 sector of Hamiltonian symmetry SU{2) ® U{1) the 
electrons read 



^Fc ' '2 ^Fs ' ' 2 ^Fc ' '2 Ms > i 2 

and the holes read 



c±fc^t,T = « J±) , _ 1 s -i ' ^ife^i.i = K^«!(±) , _ 1 «!(±) , _ 1 • (5-14) 

For the /' = 1 sector of Hamiltonian symmetry SU{2) the electrons read 



^■t,T ~ '^!{±) ^ iK'^'^!(±) , 1 ' ^ifeFia ~ ^c'^'ctji) „ iCtji) , 1 ) (5.15) 

^Fc ' '2 ys !'')2 ^Fc ' '2 ^Fs ' '2 

and the holes read 



^Fc ' ' 2 "Fs ' '2 "Fc ' ' 2 ' '2 

In the sectors of Hamiltonian symmetry C/(l) (8) SU{2) we have that 

.i"; = = = = 7r[^ - i-] . (5.17) 

In the case of the / = —1 sector of Hamiltonian symmetry ^7(1) ®SU{2) the up-spin 
electron and hole read 



' ' ^Fc 2 <?s >*>-2" 9fc ''-'"2 5s '■''2 

and the down-spin electron and hole read 



•I- ^Fc ' ' 2 ys '12 ^Fr - ' 2 ^■^ ' ' 2 

For the / = 1 sector of Hamiltonian symmetry U (1) (2) the up-spin electron 
and hole read 
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4fcfT-T = , i«I(±) , 1 ' c±fc^T-T = K V*) e i«I(±) , 1 ' (5-20) 

and the down-spin electron and hole read 



cWxa = K^«!(±) , i«j±) , 1 ' cifc^i.i = ^-^%-(±',c i«!(±) , ^ • (5-21) 

^Fc ' '2 "■^ ' '2 ^-Fc ' '2 ys 2 

Finally, for the 5*0(4) initial ground state both Eq. (49) and the following 
equation 

= -gfc^ = = = a\ - ^] , (5.22) 

hold true and we find for the electrons 



If. 



(±) 



If 



(±) 



1 5 
" 2 



a 



(±) 



'2 ''Fs '■^ 



(5.23) 



and for the holes 



c±fc^T,T = j±) , _ 1 « j±) , 1 , Cifc^i.i = j±) , _ 1 « j±) , _ 1 • (5-24) 

Hpc 2 ^Fs ' '2 ^Fc ' ' 2 ^Fs ' ' 2 

Equations ( p.23D -( p.24D reveal that removing or adding electrons from the 5*0(4) 
ground state always involves creation of pseudoholes. Furthermore, while in the case 
of the (Z, V) sectors the initial and final ground states belong in general to the same 
sector, in the case of the SO (4) ground state each of the four possible transitions 
associated with adding one up-spin or one down-spin quasiparticle or quasihole leads 
to four ground states belonging to a different (/, /') sector. If the initial ground state 
belongs to the (/') SU{2) ®U{1) sector [or to the (/) U{1) ® SU{2) sector] then 
two of the final ground states belong to the (1, /') [or to the (/, 1)] sector and the 
remaining two ground states to the (—1, V) [or to the (/, —1)] sector. 



5.3 Low-energy excitations and symmetry 

Given a ground state with electron numbers {N„, A^-o-), we find in this section that 
the set of all pseudoholes of different type which constitute, in pairs, the ] and [ 
electrons and ] and J. holes associated with the ground-state transitions {N^, N^) —> 
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{N^ ± 1, A^j.) and {N-^, Ni) — > (A^|, Ni±l) transform as tlie symmetry group of the 
Hamiltonian (|4.19| ) in the corresponding sector of parameter space. 

Equations (|3.7|) tell us that the values of S"^ and S*^ are fully determined by 
the number of c, P pseudoholes whereas the number of s, [3 pseudoholes determines 
the values of S'^ and 5*. In addition, note that the quasiparticle and quasihole 
operators (^-(pID, (|5l^)-(|5l6|), (^l8D-(^:2lD, and (^^-(^2^ involve always 



a change in the number of c and s pseudoholes by one. Moreover, when acting on the 
suitable ground state these operators change the values of S''^ and by ±1/2 and 
■^sgn{Sl)l/2, respectively, and the values of and SI by ±1/2 and ±s(7n(S'f )l/2, 
respectively. (The corresponding changes in the pseudohole numbers and in the 
values of S'^, S';, S^, and are shown in Tables ^]2| and |5]^.) The analysis of 
the changes in the pseudohole numbers could lead to the conclusion that the c, ±| 
pseudoholes have quantum numbers {S'^ = 1/2; S'^ = 0; = ±1/2; 5*^ = 0) and 
that the s, =f| pseudoholes have quantum numbers {S'^ = 0; S** = 1/2; S'f = 0; S'f = 
±1/2). If this was true the c and corresponding /? pseudohole quantum numbers 
could be identified with S'^ and S'f for S""^ = 1/2 and S'; = ±1/2, respectively, and 
the s and corresponding P pseudohole quantum numbers could be identified with 
S' and 5| for S' = 1/2 and 5f = ±1/2, respectively. (Note that 5" and do 
not refer to the whole many-electron sustem but to a single electron) However, this 
is not in general true. This holds true in the particular case of zero-momentum 
number operators. On the other hand, the above identities are also true for finite- 
momentum operators for c, ±| in the limit of zero chemical potential and for s, ±| 
in the limit of zero magnetic field, as we find below. 

In order to confirm that the above equivalences are not in general true for finite- 
momentum fluctuations, we consider the a-pseudohole fluctuation operator [0, [13 



Pa{k) = -J2 (^(^ia,f^(^g+k,a,(3 , (5.25) 

and the S';- (charge) and Sf.- (spin) fluctuation operators 

(5.26) 



k'a 



_ t 



and 



PSli.k) = }2 (^4'+k,aCk',a , (5.27) 
k',cr 

respectively. From Eqs. (|3.7|) and ( p.26|) we find 



86 



5 Excitations and Finite-Size 



pM = P5s(0) =K-N^-N^, (5.28) 

and 

p,{0)=pss{0)=N^-N^. (5.29) 

We then conclude that at zero momentum the above equivalences hold true. The 
electron numbers and Ni are good quantum numbers of the many-electron sys- 
tem. Since the exact Hamiltonian eigenstates are simple Slater determinants of 
a, /3-pseudohole levels, the numbers of a, P pseudoholes are thus required to be also 
good quantum numbers. They are such that Eqs. (|5.28|) and ( |5.29| ) are obeyed. 

On the other hand, the conservation of electron and pseudohole numbers does 
not require the finite-momentum c and s fluctuations being bare finite-momentum 
charge and spin fluctuations, respectively. For simplicity, we consider the smallest 
momentum values, k ± jf-. We emphasize that for k = acting the oper- 

ator Pa{k) onto a ground state of general form ( p.43| ) generates a sinlge-pair a- 
pseudoparticle-pseudohole excitation where the a, (3 pseudohole at g = g^'* moves 
to g = q^^^. If in c, (3 the color c was eta spin and (3 = and in the 
color s was spin and [3 = S^, we should have that pc{±^) = ps^{±jf-) and 
Ps{±^) = p5|(±-^), respectively. However, the results of the works ||12|, |13| show 
that this is not true for the sectors of Hamiltonian symmetry U{1) U{1). Al- 
tough the pseudohole summations of Eqs. ( p.7|) and (|3.26|) g ive •) and 



this does not require each c pseudohole having r] spin 1/2 and spin and each s 
pseudohole having t] spin and spin 1/2. Also, the fact that the quasiparticle or 
quasihole of Eqs. O - ([5TT|) , (50) - (§J§, ( CT) - (58), and - ([Ol has 



S" = 1/2; S' = 1/2; S^^ = sgn{S^^)l/2; = sgn{S'^)l/2 does not tell how these val- 
ues are destributed by the corresponding c pseudohole, s pseudohole, and topological 
momentum shifts. 



The studies of the works |]T2|, |T3[ reveal that for finite values of the chemical 
potential and magnetic field there is a c and s separation of the low-energy and 
small-momentum excitations but that the orthogonal modes c and s are not in 
general charge and spin, respectively |jl2|, . On the other hand, these studies have 
found that in the limit of zero chemical potential the finite-momentum c fluctuations 
become real charge excitations and that in the limit of zero magnetic field the finite- 
momentum s fluctuations become real spin excitations. In the latter limit the c 
excitations are also real charge excitations and the c and s low-energy separation 
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becomes the usual charge and spin separation [|, 0, ||. 

It follows that in the case of the SO (4) canonical ensemble the set of pseudoholes 
involved in the description of the two electron and two hole operators ( |5.23|) -( p.24|) , 
which are the c, + |; c, — |; s, + |; and s, — | pseudoholes at the pseudo- Fermi points, 
transform in the S'^ = 1/2 and S*** = 1/2 representation of the SO (4) group. More- 
over, it can be shown from the changes in the BA quantum numbers and from the 
study of the pseudohole energies that the S'^ = 1/2 and 5''' = 1/2 elementary exci- 
tations studied Essler et al. are simple combinations of one of the ground-state 
- ground-state transitions generated by the operators ( ^.23| )-( p^ ) with a single 
pseudoparticle-pseudohole process relative to the final ground state. In addition, 
the usual half-filling holons and zero-magnetization spinous can be shown to be 
limiting cases of our pseudohole excitations. For instance, the {S'^ = 1/2; = 
0; 5"^ = 1/2; SI = 0) anti holon and (5^= = 1/2; S' = 0; S", = -1/2; 5f = 0) holon 
excitations are at lowest energy generated from the SO (4) ground state by the 
operators V^"a\±) ^ and V^a\±f ^, respectively. Also at lowest energy, the two 
{S^ = 0;S' = 1/2- S^^ = 0; S', = 1/2) and {S'' = 0; S' = 1/2; S^, = 0; 5f = -1/2) 
spinous are generated from that ground state by the operators a^(±) ^ and 

a^±, , , respectively. The full spectrum of these excitations is obtained by adding 
to these generators a suitable single pseudoparticle-pseudohole-pair operator. 

Since the 5*0(4) symmetry only allows Hamiltonian eigenstates with integer 
values of S'; + Sf., the holon - spinon pairs of Eqs. (|5.23|) -( ^^2^ cannot be separated. 
This also holds true in the general case, the electron being constituted by one c 
pseudohole, one s pseudohole, and one many-pseudohole topological momentum 
shift of large momentum, as confirmed by Eqs. ( ^.4|) -( prTTD , ( ^.13| )-( pn^ ), and (|5.18|) - 



( p.21j ). Also in this case the fact that only Hamiltonian eigenstates with integer 
values of S'; + are allowed prevents these three excitations from being separated. 

As for the SO (4) ground state, we can relate the symmetry of the Hamiltonian 
( [4.ig| ) in a given canonical ensemble by looking at the pseudohole contents of the 
corresponding two electrons and two holes of Eqs. (|5^)-( |5TTl| ), ( |5.13| )-( |5?T6| ), and 
(CT)-(P^)- For instance, Eqs. (^-(p^D show that in the (/,/') sectors of 



Hamiltonian symmetry U{1)^U{1) the two electrons and two holes involve one pair 
of the same type of pseudoholes, namely the corresponding c, | and s, | pseudoholes. 
Each of these transforms in the representation of the group U{1) and, therefore, the 
set of two pseudoholes transforms in the representation of the group U{1) ® U{1). 

In the case of the (/') sectors of Hamiltonian symmetry SU{2) ® [/(I) c is r/ spin 
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and the corresponding quantum number j3 is SI and Eqs. (|5.13| )-( |5.16| ) confirm that 
the two electrons and holes involve either one c, | pseudohole or one c, — | pseudohole 
combined with one s, |- pseudohole. The c, | and c, — | pseudoholes transform in 
the S"^ = 1/2 representation of the r^-spin SU{2) group, whereas the s, |- pseudohole 
transforms in the representation of the U{1) group. Therefore, the set of c, |; c, — |; 
and s, I pseudoholes transforms in the S'^ = 1/2 representation of the SU{2) ®U{1) 
group. 

In the case of the (/) sectors of Hamiltonian symmetry f/(l) SU{2) s is spin 
and the corresponding quantum number P is P = and Eqs. (|5.18|) - (|5.21| ) show 
that the two electrons and two holes are constituted by either one s, ^ or one s, — | 
pseudohole combined with one c, | pseudohole. The s, | and s, — | pseudoholes 
transform in the = 1/2 representation of the spin SU{2) group and the c, | 
pseudohole transforms in the representation of the U{1) group. It follows that the 
set of the c, |; s, |; and s, — | pseudoholes transforms in the = 1/2 representation 
of the U{1) ® SU{2) group. 



5.4 Finite-energy topological momentum- shift op- 
erators 

The simple form of the GGS expression (|3.38|) , GS expressions (p.43|) -( ^^45|) , and 
of the general-Hamiltonian eigenstates ( |3.18| ) introduced in Chapter ^ has a deep 
physical meaning. It reveals that in the present basis these eigenstates of the many- 
electron quantum problem are "non-interacting" states of simple Slater-determinant 
form. However, that the numbers I"''^ of the rhs of Eq. (28) can be integers or half- 
odd integers for different Hamiltonian eigenstates, makes the problem much more 
involved than a simple non-interacting case. This change in the integer or half-odd 
integer character of some of the numbers /"'^ of two states, shifts all the occupied 
pseudomomenta. As in the case of last section, transition between a GS (46) and any 
eigenstate (|3.18| ) can be separated into two types of excitations: (a) a topological GS 
- GGS transition which involves the creation or annihilation of pseudoholes and (or) 
heavy pseudoparticles as well as the occurrence of topological momentum shifts and 
(b) a Landau-liquid excitation associated with pseudoparticle - pseudohole processes 
relative to the GGS. 

The topological transitions (a) are basically superpositions of three kinds of el- 
ementary transitions: (i) GS - GS transitions involving changes in the a electron 
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numbers by one, (ii) single (3 flip processes of a, (3 pseudoholes which lead to non- 
LWS's and non-HWS's outside the BA, and (iii) creation of single a, 7 pseudopar- 
ticles at constant values of S'^, 5*^, S^, and S^. While the transitions (i) are gapless, 
the elementary excitations (ii) and (iii) require a finite amount of energy. 

The generators of the excitations (i) were studied in the previous section and in 
the work [|ll| and the ones of (ii) involve only topological momentum shifts and f3 
flips (which describe either creation of electron pairs or spin flip processes). Consider 
as an example of a transition (iii) the creation of one c, 7 pseudoparticle which is 
found to be the relevant process for the finite- frequency conductivity flSl. (We 



examine here explicitly topological momentum shifts of occupied bands only, i.e. 
these which generate momentum.) When 7 > is even (or odd) the generator reads 
V^^^(jc,7 (or Gc^^V^^), where the topological-momentum-shift operator was defined 
above and 

(-) ,22^( I) (+) 

Gc,, = [sir n n ^L-i^Uc.- (5-30) 

0=0'"^ a'+' - — (-y-l) ^ 

y 'JPcfi ^Fc,0 Na \ ' ' 

Therefore, such transitions involve one c or s topological momentum shift, the cre- 
ation of a number 7 of c, | pseudoholes and 7 of c, — | pseudoholes, and the creation 
of one c, 7 pseudoparticle at g = 0. Similar results hold for the topological momen- 
tum shifts associated with the creation of one 5,7 pseudoparticle. In next section 
we show the physical importance of this latter momentum shifts. 



5.5 Momentum shifts at work: the example of 
spin waves 

By use of a simple example, we now show that the topological momentum shifts 
contribute to the total momentum of a given excitation. We consider ID Hubbard 
model at electronic and magnetization densities equal to 1 and 0, respectively, and 
in the strong coupling regime U/At^ ^ 1. 

It was shown by Ogata et al. [|T^ that in the strong coupling regime the Bethe- 
ansatz wave function for the Hubbard model decouples into the product of the ID 
anti-ferromagnetic Heisenberg model wave function and a Slater determinant 



of free spinless fermions. At n = 1 the spinless fermions have no dynamics and we 
are left with the antiferromagnetic Heisenberg chain with magnetic coupling given 
by At^ /U . The correct description of the anti-ferromagnetic chain ( p.20| ) excitation 
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T8|. These authors 



spectrum at zero magnetic field was given by Faddeev et al. 
found that the excitation spectrum is two parametric, which leads to a continumn 
of excitations. The excitation spectrum E{qi, ^2) and the excitation momentum AP 
are given by 



Jtt 

E{qi, qi) = — [cos(gi) + cos(g2)] , AP = tt - gi - 52 • (5.31) 

Here J is the antiferromagnetic coupling and qi and q2 are the momentum of two 
spin waves (each with spin 1/2). The factor vr in AP expression refers to the ground- 
state momentum (A^^^ was choosen to be odd and the number of sites Na was choosen 
to be even). 

The same results can be obtained for the Hubbard chain in the strong-coupling 
regime. Although the spinless fermions band is full, all the fermions change their 
momenta by vr/A^a, what recovers the excitation momentum of the Heisenberg chain. 
For f//4t^ ^ 1 and electronic density n = 1, the energy band of the s, pseudopar- 
ticles is given by [|19| 

eloi<l) = —^^^os{q) (5.32) 

and the c, energy band by 

£l^{q) = -2t cos(g) + 0{4tyU) . (5.33) 

At zero magnetic field the ground state is the pseudohole vacuum, that is the SO (4) 
ground state. In this case the s,7-energy bands collapse to the point e'^^^iq) = 
and q = 0. Both the c, and the s, bands are filled in the ground state and 
^^,o{lFs,o) = £^8,7(0) = 0. Thus, creation of the s,7 heavy-pseudoparticles costs no 
energy (excitations involving creation of c, 7 heavy-pseudoparticles have an energy 
gap). That is, excitations between the ground state and LWS's II. Let us consider 



that the initial ground has both and N-^ odd. From Table |5.1| we see that the 
ground state has zero momentum. The simplest excitation one can think of is, 
under these conditions, the creation of one single s, 1 heavy pseudoparticle. From 
the numbers (|3.26|) and keeping constant Ni and N-^, we see that the total number 
of s, pseudoholes varies by 2, and that the number of c, pseudoparticles remains 
unchanged (the corresponding number of pseudoholes is zero). ( Equation ( p.lO| ) 
confirms that the s, 1 pseudo-Brillouin zone is reduced to the single q = point.) 



Equations (|3.12|) and (|3.13|) imply that all c, pseudoparticles (spinless fermions 
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in this limit) change their momenta by ir/Na. This leads to a finite momentum of 
6q = TT, and the energy and momentum excitation spectra are 

E{qi,q2,q3) = ^sAqs) + ^[cos{qi) + cos{q2)] = — -[cos(gi) + cos(g2)] (5.34) 
and 



AP = TT + - gi - g2 = TT - gi - , (5.35) 

respectively. It is obvious that Eqs. ( p.34| ) and ( |5.35| ) recover Faddeev's results ( p.31| ) 
for the Heisenberg chain and that our s, pseudoparticles are Faddev's spin-1/2 spin 
waves, in agreement with the results of Sec. 

The same excitation spectra, ( |5.34| ) and ( |5.35| ), are obtained if we take the limit 
H = 0^ and make a triplet excitation (fiip one electron spin). In this limit and for 
the same parity of the electronic numbers N^j we have just used, the ground state 
has two holes in the s, band, as Eq. ( |3.10| ) implies. The ground-state momentum 
is zero. Due to the electronic spin-fiip excitation, the boundary conditions ( |3.12| ) 



and (|3.13|) adjust in such a way that the final state still has two holes in the s, 



band. The number of s, pseudoparticles changes by one, and a momentum shift 
of Tc/Na occurs in the c, band, what leads, at half- filling, to a finite momentum of 
TT. We then recover the result that at zero magnetic field the triplet and the singlet 
excitations are degenerated. 

Without the use of the Bethe-ansatz equations, Dias and Lopes dos Santos found 
equivalent results For U = oo they confirmed that all the spin configurations 



are degenerated. This corresponds, in the pseudoparticle formalism, to £^°^>o(q') = 
0. Despite this massive degeneracy, changes in the spin configurations lead to a 
global momentum shift of ir/Na in the c, pseudoparticle band. This effect was 
represented by Dias and Lopes dos Santos as a fictitious magnetic fiux through 
the spinless-fermion ring. This simple representation shows that the change in the 
pseudomomenta is such that the system minimizes its energy when both the spin 
configuration changes and one electron is added to or removed from the system. 
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Chapter 6 

Zero- Temperature Transport 



6.1 Introduction 

The transport properties of strongly correlated electron models for low- dimensional 
conductors has been a subject of experimental and theoretical interest for over 
twenty years. Low-dimensional conductors show large deviations in their transport 
properties from the usual single-particle description. This suggests that electronic 
correlations might play an important role in these systems 0, |, |], |[, even if 
the on-site Coulomb repulsion is small [^]. 

In this chapter we use the generahzed pseudoparticle theory to study the charge 
and spin currents of the Hubbard chain at finite energy. For that we solve the 
BA equations with a twist angle for all electronic densities and magnetizations. 
We express the charge and the spin-diffusion currents in terms of the elementary 
currents of the charge and spin carriers. It is shown that the latter carriers are the 
a, 7 pseudoparticles of the pseudoparticle-perturbation theory (PPT) presented in 
Chapter |^ and introduced in the work 0. We evaluate their couplings to charge 
and spin and define the charge and spin pseudoparticle transport masses. The ratios 
of these masses over the corresponding static mass provide important information 
on the role of electronic correlations in the transport of charge and spin in the 
ID quantum liquid. Furthermore, we find that the transport of charge and spin 
can be described by means of pseudoparticle kinetic equations. Our results are 
a generalization to finite energies of the low-energy results on transport of charge 
and spin results obtained in earlier works [§, |^. This is possible by means of the 
generalized pseudoparticle representation introduced in Chapters || and 

Recently, the pseudoparticle description of all BA Hamiltonian eigenstates in- 
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troduced in Chapter ^ has allowed the evaluation of analytical expressions for 
correlation functions at finite energy |1^. From these results one can obtain ex- 
pressions for the absorption-band edges of the frequency-dependent electronic con- 



ductivity, (j{uj) [115 



6.2 Review of fundamental concepts in pseudopar- 
ticle theory 

Although the pseudoparticle description, introduced in Chapter |, refers to all 
Hamiltonian eigenstates 0, in this Chapter we restrict our study to the Hilbert 
subspace involved in the zero-temperature charge and spin frequency-dependent 
conductivity |jl5[. This is spanned by all the Hamiltonian eigenstates contained 
in the states j'^\GS), where \GS) denotes the ground state and j'' are the charge 
(C = p) and spin {( = a^) current operators [given by Eqs. ( |6.15| ) and ( |6.16| below. 



respectively]. Since these current operators commute with the six generators of the 



?7-spin and spin algebras [jT6[ (Chapter pTT|) , our Hilbert subspace is in the present 



parameter space spanned only by the lowest-weight states (LWS's) of these algebras 
[IT?! (Chapter This refers to the Hilbert subspace directly described by the BA 
solution [^]. [Therefore, following the studies and notations of Chapter || 0, we 
can use the a,7 = pseudoparticles instead of the pseudoholes (with = ±|) 
required for the description of the non-LWS's outside the BA solution.] 

The c, and s, pseudoparticle branches have been called in previous low-energy 
studies c and s pseudoparticle branches, respectively [||, ||. They were shown to 
describe the low-energy excitations of the Hubbard chain and to determine the low- 
energy behavior of its charge and spin transport properties p, ^. (In the limit 
of low energy the heavy-pseudoparticle branches are empty.) On the other hand, 
description of the LWS's of the model that have a finite-energy gap, Uq, relatively 
to the ground state involves the heavy pseudoparticle branches c, 7 > and s, 7 > 
(Chapter I) [0. 

A very useful concept in this theory is that of generalized ground state (GGS). 
In Chapter |^ it was defined as the Hamiltonian eigenstate(s) of lowest energy in the 
Hilbert subspace associated with a given sub-canonical ensemble. The concept of 
sub-canonical ensemble follows from the conservation laws of the a, 7 pseudoparticle 
numbers, Na,-y Each Hamiltonian eigenstate has constant values for these numbers 
and a sub-canonical ensemble refers to a given choice of constant Na,-y numbers. 
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On the other hand, the concept of GGS has been recently extended p!3| , as 



follows (i) filled a, 7 pseudoparticle seas with compact occupations around q = 0, i.e. 

?Fa,7,+i Q ^ ?Fa,7,+i' wherc the pseudo- Fermi points are given by 5'ira',7,+i — 
_j_7r^j, _|_ (9(^]^/]Y^)^ a.nd (ii) filled 0,7 pseudoparticle seas with compact occupations 



for <q< 

'?Fa,7,-i ^'^'i ^Fa,^,-! ^ Q ^ la,] ^ whcrc the pseudo- Fermi points 
±) ^ I r _ 



are given by q^pa^-^^^i = ±[5'a,7 — ""^"'^ ] + ©(l/Ai'a). Recent developments on the 



problem [0, |T^, |I5|, have revealed that the creation of one a, 7 pseudoparticle 
from the ground state involves, to leading order, a number 27 of electrons. Since 
the currents are two-electron operators, it follows that the creation of single a, 1 
pseudoparticles from the ground state are the most important contributions to the 
transport of charge (a = c) and spin [a = s) at finite energies. On the other 
hand, since the states j''\GS) which define our Hilbert space have zero momentum 
(relatively to the ground state) and the creation from the ground state of single a, 1 
pseudoparticles of type (ii) is a finite-momentum excitation, for simplicity in this 
chapter we restrict our study to GGS's of type (i). Therefore, in order to simplify 
our notation we denote the pseudo- Fermi points qp^a^^^+i simply by qp'a,'y These are 
given by qp'a,^ = ^qFa,-y + 0(1/ No) where the pseudo- Fermi momentum (Chapter 
I) [0 

qpa^-r - — , 10- J- j 

appears in several expressions below. Note, however, that the generalization of our 
results to GGS's of type (ii) is straightfoward. We emphasize that in Chapter |^ the 
definition of GGS refers to the above choice (i) for the c, and s, 7 branches and to 
the choice (ii) for the c, 7 branch with 7 > 0. Therefore, in the case of the c, 7 > 
pseudoparticles our GGS choice differs from the choice of that chapter. 

The ground state is a special case of a GGS where there is no 0,7 > heavy- 
pseudoparticle occupancy and the pseudo-Fermi points (|6.1|) are of the form 



qpcfi = '^kp ; qpsfi = kpi ; gFa,7 = 7 > . (6.2) 

Since the conservation of the electron numbers imposes the following sum rules 
on the numbers Na,'y (Chapter |^) 

^i = E7iVc,7 + E(l + 7)iVs,7, (6-3) 

7>0 7 

and 
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iV = iV,,o + 2 5^7iVc,7, (6.4) 

7>0 

the creation of heavy pseudoparticles from the ground state at constant electron 
numbers requires the annihilation of a, pseudoparticles. It follows from Eqs. ( p.3|) 
and (|6.4|) that the changes AA^^q, o associated with a corresponding ground-state - 
GGS transition read 



s,0 



E7iVc,7-E(l + 7)iV. 

7>0 7>0 



s,7 1 



(6.5) 



and 



AiV,,o = -2^7iVc,7. (6.6) 

7>0 

For instance, the creation of one c, 7 heavy pseudoparticle from the ground state 
requires the annihilation of a number 27 of c, pseudoparticles and of a number 7 
of s,0 pseudoparticles, whereas the creation of one s,7 pseudoparticle involves the 
annihilation of a number 1 + 7 of s, pseudoparticles and conserves the number of 
c, pseudoparticles. 



6.3 Charge and spin currents: solution of the BA 
Equations 

Within linear response theory the charge and spin currents of the ID Hubbard model 



can be computed by performing a spin- dependent Peierls-phase substitution [|18|, |T9 
in the hopping integral of Hamiltonian ( |2.1(J| ), t f^±i<t>a/Na^ leading to Hamiltonian 

It has been possible to solve the Hamiltonian ( p.lO|) with the additional hop- 
pmg phase e^^'^''/^" by means of the coordinate BA both with twisted and toroidal 



boundary conditions, both approaches giving essentially the same results p!9| , |20 
One obtains the energy spectrum of the model parameterized by a set a numbers 
{kj,As} which are solution of the BA interaction equations given by 



N 



sm{kj) -As + iU/A 
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and 



sin(A;,) -As- iU/4 ~ 



A/3 - A5 - iUl2 ' 



(5 



(6.8) 



These equations reduce to Eqs. ( |2.36| ) and ( |2.37| ) with 



However, previous studies of the 0o- 7^ problem |19, 20 1 have only considered 



the real BA rapidities of Eqs. ( |6.7|) and ( |6.8|) which refer to low energy. Here we 
follow the same steps as Takahashi |^ for the 00- = Eqs. ( |6.7| ) and ( |6.8| ) and 
consider both real and complex rapidities. We then arrive to the following 0^ 7^ 
equations which refer to the real part of these rapidities 



Ns. 



k,Na 



27r/| + 0T - I] II 2 tan 

7 j'=l 



I sm{kj)/u — Rs,'y,j' 



(7 + 1) 



EE2tan 

7>0 j'=l 



_^ / sin(A;j)/'u — Rc,-y,j' 
V 7 



(6.9) 



2NaRe sin-i([i?e,7j + ^7]^) = 27rlf"' + 7(0| + 0J 

- E2tan-f^^^fel/^i^^) 
i'=i V 7 / 

+ E E ®7,7'(-^c,7,j - -Rc,y,j') > (6.10) 
7'>0 j'=l 



and 



, /i?,7,- -sin(A;')/M\ , 
E 2tan-i (1 + ^) J + (1 + 7)(0T - 0i) = 

= 2vr/;'^ + E E ©7+i,7'+i(^s,7J - ^s,7'y) • (6-11) 
7' i'=i 

The functions 67,7/ (x) [and 67+i_7/+i(a;)] of Eqs. (|6.9|), (|6.10|) , and (|6.11|) are defined 
by Eq. (|4.8|) . The following definitions for the real part of the rapidities, A'^'^^/u = 
Rs,'y,j (with n + 1 = 7 and a = j), A'J'/u = Rc,^j (with n = 7 and a = j), and 
7 = 1, 2, . . . , 00 for the Nc^^ sums and 7 = 0, 1, 2, . . . , cxd for the A^s,7 sums, allows 
us to recover Takahashi's formulae for = [0]. Here and often below we use the 
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notation c = c, 0, which allows the c, 7 sums to run over 1,2,3,..., 00. Whether we 
are using this notation or the previous one will be obvious from the context. 
The important numbers H, If^, and Ij''^ which appear when going from Eqs. (|6 



and (6^) to Eqs. (|6.9| ), ( |6.10 ), and ( |6.11 ) are the quantum numbers which describe 



the Hamiltonian eigenstates. As we have seen in Chapter ^ it is convenient to de- 
scribe the eigenstates of the model in terms of pseudomomentum {q"'^ = 2ttI"'^ /Na} 
distributions, where = Ij. The energy and momentum eigenvalues are given by 



j=l 7>0 j=l 

+ NaiU/A-^) + N{fi-U/2)-^oHiN^-N^), (6.12) 



and 



Ns 



E(^; + f) + EEar-(7 + i) 



7 j=i 



7>0 j=l 



+71 J2 , 

7>0 

(6.13) 



respectively. Eqs. ( |6.12D and (|6.13|) are the eigenvalues of Eqs. ([4.19|) and ( p.32| ). 
We emphasize that the rapidity dependence on (p^ is defined by Eqs. (|6.9|)-( |6rTlD 
and determines the energy-functional ( |6.12| ) dependence on 0^. The corresponding 
0O- = expressions recover the rapidity equations (|6.9| ), ( |6.10| ), and (|6.11| ). 

In the limit of a large system (A^^^ 00, N/Na fixed), and using the same 
methods as in Chapters ^ and ^, we can develop a generalization of the low-energy 
pseudoparticle-Landau-liquid description of the Hubbard model. This generalization 
refers to energies just above the set of energies c^o given by 



2/i^7iVc,7 + 2/ioi^E(l 

7>0 7>0 



E <7(o)^. 



0,7 ) 



(6.14) 



o,7>0 



where the Hamiltonian ( |4.40| ) truncated up to second order describes the quantum- 
problem at {uj — ujq) low energies. Note that the choice Uq = 0, which refers to 
A^a,7 = for 7 > 0, recovers the usual low-energy theory [|, On the other hand, 
when uq > 0, in addition to finite occupancy of the usual c, = c and s,0 = s 
pseudoparticle bands, there is finite occupancy for some of the branches of the c, 7 
and s, 7 heavy-pseudoparticles M . 
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In the above thermodynamic hmit the rapidity real parts kj = kj{qj), Rs,-y,j = 
Rs^^j{qj), and Rcyj = -Rc,7,j(?j) prohferate on the real axis. As we have done in 
Chapters |^ and Eqs. ( |6.9|) , (|6.10|) , and (|6.11|) can be rewritten as integral equa- 



tions with an explicit dependence on the pseudomomentum distribution functions 
Na^:y{q). These are Eqs. ( |A.7|) , ( |A.8|) , and ([A.9|) of Appendix ^ which refer to the 



case 0O- 7^ 0. In that Appendix we derive ground-state normal-ordered expressions 
for the rapidities and charge and spin currents. 

The combination of Eqs. (|6.12|) , (|A.7|) , (|A.8| ), and (|A.9|) allows the evaluation of 



several interesting transport quantities. This includes the charge and spin currents 
and the charge and spin pseudoparticle transport masses. The charge and spin 
current operators (with C, = p for charge, and C = ^2 fo^' spin) are for the ID 
Hubbard model given by |[ 

Na 

jP = -ezt ^ ^(c]^Cj+i^ - c]+i^Cj^) , (6.15) 

a j = l 

and 

Na 

r = -(l/2)zt5:^a(49+i.-c]+,,c,.). (6.16) 

cr j = l 

Importantly, the discrete nature of the model implies that the commutators of the 
Hamiltonian ( |2.10 ) and of the current operators j^, Eqs. ( |6.15| ) and ( |6.16 ), are non 



zero. It follows that the BA wave function does not diagonalizes simultaneously 
the Hamiltonian ( p.lO|) and the current operators (|6.15|) and ( |6.16D . Since the BA 
solution alone only provides the diagonal part in the Hamiltonian-eigenstate basis of 
the physical operators [0 , we can only evaluate expressions for the diagonal part of 
the currents which provide the mean values of the charge and spin currents. These 
refer to all LWS's and are important quantities for they allow us to compute the 
transport masses of the charge and spin carriers of the system. In addition, our 
formalism defines the charge and spin carriers. These are found to be the c (i.e., 
c, 0) and c, 7 pseudoparticles for charge, and the c and s, 7 pseudoparticles for spin. 
This follows from Eqs. ( |A.7| )- (|Al9| ) and also from the Boltzmann transport analysis 
of Sec.|6.5|. 

We emphasize that combining the generalized pseudoparticle representation of 
Chapter |^ with a low-energy {oj — ujq) conformal-field theory [jl3|, leads to 



finite-energy current - current correlation function expressions which are determined 
by the non-diagonal terms (in the Hamiltonian-eigenstate basis) of the current opera- 
tors. This is a generalization of the low-energy correlation-function studies of several 



102 



6 Zero-Temperature Transport 



authors |T0|, |TT], |2^, However, these expressions cannot be derived within the 
BA solution alone. Therefore, these studies go beyond the scope of the present the- 
sis and here we consider the diagonal part of the charge and spin current operators 
only. 

The mean value of the current operator evaluated for the LWS \m) is given 

by 



(m| y\m) = — 



(6.17) 



(/)=0 



where Em is the energy of the Hamiltonian eigenstate \m) and |T 



= 0t = -0i, C = ^z. (6.18) 

In our basis the LWS's are simply obtained by considering all the possible occu- 
pation distributions of the pseudomomenta qj = 2t:I'^'^ /Na- Therefore, it is conve- 
nient to describe the matrix elements (m| j^|m) in terms of the pseudomomentum 
occupation distributions Na,-y{q)- This leads to a functional form for the current 
mean values. The computation of {m\j''\m) involves the expansion of Eqs. ( |6.12D 



and ([A.7|) - ([A.9|) up to first order in the fiux (p. Writing Eq. (|6.12| ) in the limit of 
Na ^ oo, expanding it up to first order in the fiux 0, and using Eq. (|6.17|) we obtain 



(m| jV) = -^twz I dqN,{q)K'^{q)sm{K{q)) 

-qc 



2tt 



+ E «^ ^,iV„(,)fle RtM , (6.19) 

7>0 -'-I-'-' ^1 - U^[Rca{<l) - ^l] 

where the important functions W^i^q) (with W = K, Rg^-y, and Rc,y) are the deriva- 
tives of the rapidity functions defined by Eqs. ([A. 7]) - (|A.9D in order to the fiux (p at 



= 0. They obey a set of integral equations obtained from differentiation of Eqs. 



It is convenient to write (m| p |m) in normal order relatively to the ground state. 
To achieve this goal we expand all the rapidities W{q) and the functions W'^{q) as 

W{q) = W°{q)+W\q) + ... , (6.20) 
W^{q) = W°''l'{q) + W^'^iq) + . . . , (6.21) 
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respectively. In these equations the functions W^{q) and W^''^{q) are both referred 
to the ground state, and the functions W^{q) and W^'^^i^q) are first-order functionals 
of the deviations 6Na^y{q). In Appendix ^ we show that the above expansions lead 
to a ground-state normal-ordered representation. To first order in the deviations, 
the normal-ordered expression for the matrix element ( p.l9| ) simply reads 



(m| j^|m) =J2J2[ ^ dq6Nc,,j{q)ji^^{q) 



where the elementary-current spectrum ji^{q) is given by 



(6.22) 



(6.23) 



Here 



Ztt 



a'r('[QjJQFa',-y') , 



j=±l 



(6.24) 



and Ci are the coupling constants of the pseudoparticles to charge and spin given 

by 



^a,7 ^a,c^j,0 ~l~ ^a,7 ; 



(6.25) 



where 



/C^,^ = <5„,,27 ; /C^;^ = -Sa,s2{l + 7) • (6.26) 

As in a Fermi liquid p4|, ^ , the expressions of the elementary currents (|6.23|) in- 
volve the velocities Va^^yiq) and the interactions [or /-functions] fa,r,a',Y{lj'l')- How- 
ever, the pseudoparticle coupling contants to charge and spin, Eqs. ( |6.25| ) and ( |6.26|) , 
are very different from the corresponding couplings of the Fermi-liquid quasiparti- 
cles. We emphasize that at low energy Eq. ( |6.22| ) recovers the expression already 
obtained in previous works @, Q which only contains the c = c, and s = s,0 el- 
ementary currents. The coupling constants (|6.25|) - (|6.26|) play an important role in 
the description of charge and spin transport and are a generalization for 7 > of 
the couplings introduced previously p], |^. They define the a, 7 pseudoparticles as 
charge and spin carriers. We emphasize that when = the corresponding 0,7 
pseudoparticles do not couple to ( (i.e. charge or spin). Therefore, for 7 > the 
c, 7 and s, 7 pseudoparticles do not couple to spin and charge, respectively. This 
is related to the charge and spin separation of one-dimensional quantum liquids. 
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Importantly, when C^^ = the a, 7 pseudoparticle - pseudohole processes do not 
contribute to the ( correlation functions. 

In contrast to the general current expression (|6.19| ), expression ( |6.22D is only 
valid for Hamiltonian eigenstates which differ from the ground-state pseudoparticle 
occupancy by a small density of psudoparticles. This is because in expression ( |6.22|) 
we are only considering the first-order deviation term. 

The velocity term of the current-spectrum expression (|6.23| ) is what we would 
expect for a non-interacting gas of pseudoparticles and the second term takes ac- 
count for the dragging effect on a single pseudoparticle due to its interactions with 
the remaining pseudoparticles. (This is similar to the Fermi-liquid quasiparticle el- 
ementary currents [|^, |2^.) We remind that Eq. (|6.23| ) is valid for finite energies 
LO just above the energy c^o corresponding to the suitable set of finite Na^^y num- 
bers. These numbers characterize the state |m). Therefore, the sum over 7 is in 
Eq. ( |6.23|) restricted to the a, 7 bands that have non-zero occupancy of pseudopar- 
ticles, as is imposed by the presence of the step function. Within the PPT, the 
deviation second-order pseudoparticle energy expansion corresponds to the devia- 
tion first-order current expansion ( |6.22| ) which refers to small positive values {u — ujq) 
of the excitation energy. In contrast to Fermi liquid theory, our PPT is valid for 
finite energies [just above the energy values uq, Eq. ( |6.14|) j because (i) there is only 
forward scattering among the pseudoparticles at all energy scales and (ii) at low 
(a; — u;o) energy values only two-pseudoparticle forward-scattering interactions are 



relevant. (In a Fermi liquid this is only true for luq = and u ^ [0, |25[.) We 



emphasize that the current expression (|6.19|) includes all orders of scattering and, 
therefore, applies to all energies without restrictions. 



6.4 Pseudoparticle static and transport masses 



In previous studies P, || the charge and spin transport masses of the c, and s, 
pseudoparticles were defined and were shown to play an important role in the trans- 
port of charge and spin. For instance, they were shown to fully determine the charge 
and spin stiffnesses p, |^, |19|, ^ Here we generalize the mass definitions 



of these studies to 7 > and define the charge and spin transport masses, m£ , as 
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where = ji^y{qFa,'y) ■ These masses contain important physical information. As 
in a Fermi hquid p4|, E5l, the ratio 



a, 7 



^a,7/''^a,7 ) 



(6.28) 



of the transport mass, m^ ,^, over the static mass, wi* provides a measure of the 
electronic correlation importance in transport. Similarly to the 7 = case p, ^, 
the latter is in general defined as 



m 



«,7 



QFa,- 



(6.29) 



In Appendix |B| we define the mass ( p.29| ) in terms of suitable functions and find 
some limiting expressions. 
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Figure 6.1: The ratio m^^/m*^ as function of U at electronic density n = 0.7 and 
for values of the magnetic field h = H/Hc = 0.1 (full line), h = 0.5 (dashed line) 
and h = 0.9 (dashed-dotted line). For other electronic densities, the plots follow the 
same trends as for n = 0.7. 

It can be shown from the transport- and static-mass expressions that the ratio 
Tn^ /m^ involves the Landau parameters 



a,7;a' ,7' 



a,7;Q ,7 



',7'(9-Fo,7) ; 



0,1 



(6.30) 



with F^^^.a',yil) given by Eq. (|6.24|) . These parameters can be written as follows 
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Fa,r,a',Y - -^a,a'S^,YVa,^ + ^ ^(^a",7")^a",7" [C",7";a,7C",7";a',7'] ' (6.31) 



a", 7" 

where the quantities Ca,r,a','y' given by 



(9Fa,7, ^Q'Fa',7') • (6.32) 

l=±l 



We find for the ratios rn'^^/m'^^ the following expressions 



m 



a,0 



7 = 0, (6.33) 



and 



m' 



C 

a, 7 



"^a,7 Ca,7(Ca,7 + X^a' ^a',0^a,7;a',o) 



7>0. (6.34) 



In the Table ^TT] analytical limiting values for the mass ratios of form (|6.28|) are 
listed. Obviously, since for 7 > the c, 7 and s, 7 pseudoparticles do not couple to 
spin and charge, respectively, the ratios m^^^/m* ,^ and m^,^/m* ,^ are infinite. 



H ^ H ^0 n 



^cal^la 27(27-r,^_i) 27(27+i<i,^.,,o) 

2(7+1) (27+2-r,^) 2{7+l){27+2+2e_^.^.o) 2{7+l){27+2+2Ci_^.. ^ q) 



Table 6.1: Mass ratios in several limits of physical interest. The function 77^ is defined 
as Tj^ = 2/(7r) tan~-'^[(sin(?27r))/(M[7 + 1])]. The equations for the static masses rn^ ,^ 
are given in Appendix |B[ In the case H —>■ 0, simple expressions for the parameters 
^a,7;a',0' ^'i- ( |6-32| ), cau be obtained from the results of Appendix |^. The ratios 
m^^^/m*^ and mP^/m*^ are infinite. The dependence on U and n of the parameter 
^0 has been studied by Frahm and Korepin [E 



As was referred previously, it can be shown from the results of Chapter ^ and 
other studies |0, [l^, that the creation of one a, 7 pseudoparticle from the 
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Figure 6.2: The ratio m'^i/m* ^ as function of the electronic density n and for values 
of the magnetic field h = 0.1, h = 0.3, h = 0.5, h = 0.7, and h = 0.9. The onsite 
Coulomb interaction is (a) f/ = 1 and (b) U = 5. 



ground state is a finite-energy excitation which, to leading order, involves a number 
27 of electrons. Therefore, and since the current operators are of two-electron char- 
acter and couple to charge and spin according to the values of the constants ( |6.25 ) 
and ( |6.26| ), at finite energies the c, 1 and s, 1 heavy pseudoparticles play the major 
role in charge and spin transport, respectively. On the other hand, the a,7 > 1 
heavy pseudoparticles have a minor role in charge and spin transport. It follows 
that in the present section we restrict our study to the ratios (|6.28|) involving 7=1 
heavy pseudoparticles. We consider the ratios Tn^c,il^*c,i ^"^^ "^sfi/^^s.i- (Note that 
^c,i/"^c,i = ^s,i/^*s,i = ^0 We also consider the case of the c, charge-mass ratio 
which is closely related to the charge stiffness studied in detail by other authors 
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Figure 6.3: The ratio m'^ i/m* ^ as function of the the magnetic field h and for values 
of the electronic density n = 0.1, n = 0.3, n = 0.5, n = 0.7, and n = 0.9. The onsite 
Coulomb interaction is (a) [/ = 1 and (b) U = 5. 



||. In Figs. |6?T1 - |6.12| these ratios are plotted as functions of the onsite repulsion U 
in units of t, electronic density n, and magnetic field h = H/Hc- Note that the ratios 
of the figures are smaller than one. Moreover, the a, 1 mass ratios never achieve 
the value 1 whereas the a, mass ratios tend to one in some limits because of the 
generalized adiabatic principle [§, |]. 



Combined analysis of Figs. ^3] and |6]^ reveals that the charge-mass ratio for 
the c, 1 pseudoparticle is, for large U, fairly independent both of the band filling 
n and magnetic field h. It is a decreasing function of U and as a function of the 
density, n, goes through a maximum for a density which is a decreasing function of 



U. Moreover, figures 3.3 and 6.4 show that this ratio is a decreasing function of the 



6.4 Pseudoparticle static and 



109 




(b) 



Figure 6.4: The ratio m'^i/m* ^ as function of the the magnetic field h and for values 
of the onsite Coulomb interaction U = 1, U = 2, U = 3, U = 5, U = 10, and U = 20. 
The electronic density is (a) n = 0.5 and (b) n = 0.9. 



magnetic field. 



In contrast, Figs. |6]^ - 16]^ reveal that the charge-mass ratio for the c, pseu- 
doparticle is an increasing function of U and of h and as a function of the density, 
n, goes through a minimum for a density which is a decreasing function of U. Note 



that from Fig. |6.5| the evolution of the c, pseudoparticles to free spinless fermions 
as U increases is clear. This is signaled by the ratio going to one as U oo. This 
behavior follows from the generalized adiabatic principle P, ^ and agrees with the 
well known decoupling of the BA wave function in free spinless fermions (in the low- 



energy sector and localized antiferromagnetic spins [^. Figures |6.?| and |6]| 
also reveal that in the fully-polarized ferromagnetic limit, h 1, the ratio goes to 
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Figure 6.5: The ratio Tn'^ca/^co function of U , for electronic density n = 0.7, and 
for values of the magnetic field h = 0.1, h = 0.3, h = 0.5, h = 0.7, and h = 0.9. For 
other electronic densities, the plots follow the same trends as for n = 0.7. 



one. This mass-ratio behavior also follows from the generalized adiabatic principle 
1^, ^ and confirms that in that limit the onsite Coulomb interations play no role in 
charge transport (they are froozen by the Pauli principle) . 

Note that in the large-?/ Figs. |6.2| - (c) and |6.6| - (c) the ratios m^^^/m*^ and 
^c,o/^*cQi respectively, are almost symmetric around the density n = 0.5. This 
implies that for large U the charge transport properties show similarities in the 
cases of vanishing densities and of densities closed to one. 



Figure |6.9| shows that the spin-mass ratio of the s, 1 pseudoparticles is a decreas- 
ing function of U but that it depends little on U for U > 6. For large U this ratio 
almost does not depend on the density n, as revealed by Fig. |6.10| - (c). Figures 
5.10| show that, in general, it is an increasing function of n but that for /i ^ 1 it 



has a maximum for an intermediate density. In figures |6.11| and |6.12| this spin-mass 
ratio is plotted as a function of h. It is a decreasing function of h. 

The transport masses are very sensitive to the effects of electronic correlations, 
as for instance to the metal-insulator transition which occurs at zero temperature 
when n —>■ 1 As a direct result of this transition, m^g ^ C)0 as n ^ 1, as was 
shown and discussed in the works p], P). Moreover, the zero-temperature charge and 
spin stiffnesses, D^, p|, ||, |T9|, ^ defined as 
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Figure 6.6: The ratio ^c.o/'^co function of the electronic density n and for values 
of the magnetic field h = 0.1, h = 0.3, h = 0.5, h = 0.7, and h = 0.9. The onsite 
Coulomb interaction is (a) f/ = 1 and (b) U = 5. 



where 



is defined for charge ( 



2 t/(0/iVj2 
p and spin 



(6.35) 



(/>=0 



cr^ in Eq. (|6.18|) , are such that 



J2alFafi/^ao ^ §]• For charge, m^Q = oo, and the latter expression reads 



27tD^ = qpcfi/^cfi ^iid is such that ^ as 1, satisfying Kohn criterion |]32 
These quantities can be computed within the pseudoparticle formalism by direct 
evaluation of Eq. ( |6.35| ). They can also be obtained by combining a pseudoparticle 
Boltzmann transport description with linear response theory P, 0. In order to con- 
firm the validity and correctness of our formalism, in Appendix we have recovered 
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(b) 



Figure 6.7: The ratio m^Q/m*^ as function of the magnetic field h and for values of 
the electronic density n = 0.1, n = 0.3, n = 0.5, n = 0.7, and n = 0.9. The onsite 
Coulomb interaction is (a) [/ = 1 and (b) U = 5. 

the charge and spin stiffness expressions (135) - (137) of the studies Q by direct 
use of Eq. ( |6.35|) . 

Equations ( |6.25| ) and ( |6.26| ) show that the Hubbard-chain charge carriers are 



the c, 7 pseudoparticles. In contrast to the zero-temperature limit where the c, 
pseudoparticles fully determine the charge stiffness, we expect that the c, 7 heavy 
pseudoparticles play an important role in the charge-transport properties at finite 



temperatures |3J, |35 



In Appendix we present simpler equations to define the pseudoparticle bands 
and phase shifts in the limit of zero magnetic field. These results show that for 
H —>■ and 7 > the bands eg^(g) collapse to a point. This is because both the 
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Figure 6.8: The ratio m'^Q/m*Q as function of the the magnetic field h and for values 
of the onsite Coulomb interaction U = 1, U = 2, U = 3, U = 5, U = 10, and U = 20. 
The electronic density is (a) n = 0.5 and (b) n = 0.9. 



bandwidth [see Eq. ( D.l )] and the momentum pseudo-Brillouin zone width (in the 
ground state) 



IT : 



k 



FT 



gc,7 = 



TT - 2kF 7 > , 
= kpi — kpi 7 > 



(6.36) 



go to zero as if — 0. This behavior is also present in the Heisenberg chain and, 
therefore, in that model the triplet and singlet exitations are degenerated at zero 



magnetic field ||3^. This also holds true for the Hubbard chain at if = and in the 
limit U ^ 4t, where the BA wave function factorizes in a spinless-fermion Slater 
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n=0.7 




Figure 6.9: The ratio rrf/i/m* ^ as function of f/, for electronic density n = 0.7, and 
for values of the magnetic field h = 0.1, h = 0.3, h = 0.5, h = 0.7, and h = 0.9. For 
other electronic densities, the plots follow the same trends as for n = 0.7. 



determinant and in the BA wave function for the ID antiferromagnetic Heisenberg 
chain. On the other hand, in the limit n — 1 the bands e^^{q) (for 7 > 0) collapse 
to a point also because both the bandwidth and the momentum pseudo-Brillouin 
zone width [see Eq. (|6.36|) ] go to zero in that limit. 



6.5 Kinetic equations for the pseudoparticles 

In the previous sections the quantum-liquid physics for energies just above the uq 
values, Eq. ( |6.14| ), was described in terms of homogeneous pseudoparticle distri- 
butions. The pseudoparticles experience only zero-momentum forward-scattering 
interactions at all energy scales. This property is absent in Fermi-liquid theory 
where it holds true only at low excitation energy when the quasiparticles are well 



defined quantum objects p^, ^ ^ This unconventional character of integrable 



models PO, HD| allows us to extend the use of the kinetic equations to energy scales 



just above the ujq energy values, Eq. ( |6.14|) , and not only to low energies 



The results presented in this section are a generalization of the kinetic-equation 
low-energy studies [^, ^. 

In the final Hilbert subspace of energy lo relative to the initial ground state 
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Figure 6.10: The ratio mj\/m*^ as function of the electronic density n and for 
values of the magnetic field h = 0.1, h = 0.3, h = 0.5, h = 0.7, and h = 0.9. The 
onsite Coulomb interaction is (a) t/ = 1 and (b) U = 5. 



the Hubbard model can be mapped onto a continuum field theory of small energy 
{u — ujq) |jT3|, 14 1 . The time coordinate t of such theory is the Fourier transform 



of the small energy {u — ujq) which corresponds to a finite energy u in the original 
Hubbard model. The validity of this approach is confirmed by the fact that it fully 
reproduces the rigorous results of Section 

Let us consider excitations described by space and time dependent pseudoparticle 
distribution functions, Na,^{q,x,t), given by 



iV„,^(g, X, t) = N^,^{q) + 5N^.,{q, x, t) 



(6.37) 
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Figure 6.11: The ratio m1\/m*^ as function of the magnetic field h and for values 
of the electronic density n = 0.1, n = 0.3, n = 0.5, n = 0.7, and n = 0.9. The onsite 
Coulomb interaction is {a) U = 1 and (b) U = 5. 

where N^^^{q) is the ground-state distribution. It follows from the PPT introduced 
Chapter || 0] and discussed in the Chapter || that the single-pseudoparticle local 
energy is given, to first order in the deviations ^A^^q, -^(g, x, t), by 

1 ■ fla'.-y' 

ea,^{q,x,t) = ea,-y{q) + — 22 dq'6Na',y{q',x,t)fc,,r,a',Y{q,q') ■ (6.38) 

Let A'^ represent the total charge, ( = p, or spin, ( = az- It follows from the 
relations ( |6.3|) and ( |6.4D involving the pseudoparticle and electron numbers that 
A''' depends linearly on the pseudoparticle deviation numbers. Thus, in the case 
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Figure 6.12: The ratio m'^y/m*i as function of the the magnetic field h and for 
values of the onsite Coulomb interaction U = 1, U = 2, U = 3, U = 5, U = 10, and 
U = 20. The electronic density is (a) n = 0.3 and (b) n = 0.5. 



of inhomogeneous excitations described by Eq. ( |6.37| ) the corresponding expectation 
value at point x and time t, (^''(x, t)), can be written as 



(6.39) 



where = —e and a'^"' = 1/2. 

In this "semi-classical" approach the response to a scalar field, t), is pro- 

portional to the conserved quantity A'^. As for low energy in the presence of 

the inhomogeneous potential the force J-''^{x, t)a,-y that acts upon the a, 7 pseudopar- 
ticle is given by J-'^.y{x,t) = —[dV''{x,t)/dx]C^^^ x a''. It follows that the deviations 
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6Na,-y{q,x,t) are determined by the solution of a system of kinetic equations (one 
equation for each occupied a, 7 branch) which reads 



dNa^^{q,x,t) ^ dNa^^{q,x,t) dec^^^{q,x,t) _ dNa^^{q,x,t) dea^^{q,x,t) 

dt dx dq dq dx 

dN^,^{q,x,t) dV'^{x,t)^ 



dq dx 



X «^ • (6.40) 



Introducing Eq. ( p.37| ) in Eq. ( 6.40| ), expanding to first order in the deviations 



SNa,j{q,x,t), and using Eq. (|6.38|) we obtain the following set of linearized kinetic 
equations 

d5Na,y{q,x,t) ^ ^ ^d5Na,^{q,x,t) 



dt dx 
d6N^,^{q,x,t) 1 dVi{x,t) ^^^ 

dq \ dx 



, ^ /•gg'.V ^ , dSNa>^y{q',X,t) / I 



«',7' 



The conservation law for {A'^{x,t)) leads in one dimension to 

d{AHx,t)) {JHx,t)) 

dt +^^ = °' ^'-^'^ 
where {A''{x,t)) is given by Eq. ( lOoD and {J'^{x,t)) is the associate current. Mul- 



tiplying Eq. ( |6.41| ) by C^.^ x a'', summing over a and 7, and integrating over q we 
find for V^lxjt) = and by comparing the result with Eq. (|6.42| ) that the current 



spectrum ji^y{q) is given by a'' times expression (|6.23|) . (This expression has been 
derived from the solution of the BA equations with a'' = 1.) 

This agreement confirms the validity of the above low-(co' — cjq) continuum- 
field theory. The unusual spectral properties associated with the zero-momentum 
forward-scattering character of the pseudoparticle interactions follow from the inte- 
grability of the Hubbard chain & ^ 



6.6 On the conductivity of the Hubbard chain 



Solvable one- dimensional many-electron models such as the Hubbard chain are often 
used as an approximation for the study of the properties of quasi-one-dimensional 
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conductors |1|, ^, ^, |[ . Although the model has been diagonalized long ago ||21|, ^ , 
the involved form of the Bethe-ansatz (BA) wave function has prevented the direct 
calculation of dynamic response functions, including the charge-charge and spin-spin 
response functions and their associate conductivity spectra. 

Information on low-energy expressions for correlation functions can be obtained 
by combining BA with conformal-field theory |]22| , On the other hand, several 
approaches using perturbation theory |18[, bosonization [^, the pseudopar- 

ticle formalism [§, ^, scaling methods ||2^, and spin- wave theory |^ have been 
used to investigate the low-energy transport properties of the model away from 
half filling and at the metal - insulator transition BTI. Some information on the 



transport properties at finite energies has also been obtained by numerical methods 



27, 46, 47 



To study the optical conductivity of the Hubbard chain we start from its Lehmann 
representation, whose real part, 7lea{uj), is given by |19|, ^ 



TZe a(uj) 



27re^ 

h 



1 



[DS{huj) + Crreg(u;)] , 



with the charge stiffness D given by 



(6.43) 



D = — 



E 



\Hm\' 

Em — Eq 



(6.44) 



and T = — 2t X^fc.o- cos(A;)c|. ^Cfc^o- being the model kinetic energy (all other quantities 
have been previously defined in the Chapter |^) . 

Following Kohn and our previous discussion, the value of charge stiffness 
D is a criterium for the metal-insulator transition. For such transition to occur, D 
must change from non-zero to zero. According to Eq. ( |6.44|) , this is possible only 
if the two terms in this equation cancel indentically. For correlated systems with 
continuos spatial symmetry, the total current operator j'' = J2iPi/^ = P/m is 
a constant of motion (m is the particle mass, pi is the momentum of the particle 
i, and P is the total momentum). It follows that [H,j^] = 0, which implies that 
j^\m) = Pm/m\m) and 



TZe a(uj) 



2ne' 



DS{uj) , 



D 



2N, 



■(0|f|0). (6.45) 
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Equation ( |6.45| ) then implies that for systems with continuos spatial symmetry a 
metal-insulator transition can not occur and the optical-conductivity is reduced to 
the uj = absorption peak. That is, the conductivity for a correlated system with 
continous spatial symmetry behaves like a independent-electron system. 

The situation is different for lattice systems. In this case it may happen that 

is no longer a constant of motion. Then {0\j'^\m) ^ for a set of Hamiltonian 
eigenstates \m) and the optical conductivity shows a non-zero regular spectrum 
((Treg(cu) 7^ 0) and a metal-insulator transition may occur. Metal-insulator transi- 
tions can be induced, for instance, by the Peierls mechanism due to the dimerization 
of the lattice (the original band is splitted into two with the formation of a gap, due 
to the electron- phonon interaction), by Anderson localization, due to electron scat- 
tering with random distributed impurities, or by the Mott mechanism, associated 
with charge localization produced by electronic correlations. The latter mechanism 
occurs in the one- dimensional Hubbard model at half-filling. 

In the pseudoparticle description of the Hubbard model, the metal-insulator 
transition is associated with the existance of a gap ujq between the c, and the c, 7 
bands. As we have remarked before, at half-filling and zero temperature the charge 
stiffness is zero and the real part of the optical conductivity, Rea{u), has a finite 
energy gap for all values of U/ 1. For this band filling the c, 0-pseudoparticle band 
is full and charge excitations can only occur between the low-energy a,0 band and 
the c, 7 heavy pseudoparticle bands. The form of the regular part of the optical 
conductivity is determined by this type of excitations. 

As we have discussed in Section to populate, for example, the lowest c, 7 
energy-band (the c, 1 band) with one c, 1 heavy pseudoparticle one has to annihilate 
two c, and one s, pseudoparticles. The energy involved in this process is exactly 
e°i(0)-e°o(2A;^)-e°o(-2A;^)-e°o(A;^l) = 2/i(n = 1), where fi{n = 1) stands for the 
chemical potential in the limit of half-filling and €^^{q) is the pseudoparticle energy 
dispersion. In addition, this excitation process is characterized by a topological 



excitation of the s pseudoparticles of momentum —kpi ensuring that the overall 
excitation has zero momentum, as required for the transitions induced by the current 
operator (Chapter |^). In general, the process of creating one c, 7 pseudoparticle has 
an minimum energy cost of 27/2(1) and a maximum cost of 27/2(1) + 87^. Thus, 
in principle, these transitions should show up in (Xreg as a set of optical bands. 
As remarked in Section |6.2| , it is possible to show that the creation of one a, 7 
pseudoparticle from the ground state involves, to leading order, a number 27 of 
electrons. Since the currents are two-electron operators, it follows that the creation 
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Figure 6.13: Optical conductivity results for 32-site Hubbard chains at half filling 
for values of U/t = 3 (top) and 4 (bottom), computed using maximum-entropy 
analytic continuation of quantum Monte Carlo data. The solid and dashed vertical 
lines indicate the theoretical predictions for the lower edges of the first two bands. 



of one a, 1 pseudoparticle from the ground state is the most important contribution 
to the transport of charge at finite energies. Moreover, for U ^ 4t the creation 
of one c, 1 pseudoparticle is equivalent to the creation of a single doubly occupied 
site. Thus, for these values of U only a single optical band should appear in a^eg- 
This is an exact result which can be confirmed by numerical studies in this regime 

mp7i,p[. 



In order to investigate the presence of more than one optical band for small and 
intermediate values of U, recent numerical studies have been performed |jl5| , by use 
of a quantum Monte Carlo technique for systems with 32 sites. This corresponds 
to a considerably larger system than in previous exact diagonalization studies p^. 



46| . The current-current correlation function has been computed in imaginary time 
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and continued to real frequency using the so called maximum-entropy method 
The resolution of this method is limited but, nevertheless, gives a semi- quantitative 



confirmation of a multi-band optical conductivity. In Fig. |6.13| , the results are shown 
for U/t = 3 and 4 where two peaks are clearly resolved. 

As indicated in the figure, the lower edges of the spectra agree with the predicted 
gaps (Chapter H), and the second peak appears above the predicted lower edge of 
the second band, as computed from the pseudoparticle theory. This theory also 
allows the evaluation of the critical exponents associated with the conductivity edges 
1|,|1|, Note that the weight of the first band is concentrated towards the lower 



end of the allowed band of width 8t. The second peak seen in these results is 
probably dominated by the second optical band, but likely contains contributions 
also from the tail of the first band, as well as from higher bands that cannot be 
resolved due to the limitations of the method. 



Appendix A 



Normal-ordered solution of the 
BA equations with a flux 



In this appendix we derive the normal-ordered BA equations required for the eval- 
uation of Eqs. and ( lOBQ . Writing W^{q) from Eq. ( |OTD as 



= (A,l) 

dq 

where L equals Lc,o, -^c,7, or L^,^ when W equals K, Rs^-y, or -Rc,7, respectively, we 
find that W^''^{q) obeys the following equahty 

Introducing the above equation in Eq. ( |6.19| ) and writing the distributions func- 
tions Na,y{q) as N^^^{q) + SNa^-yiq), we can expand J = {m\j''\m) in terms of the 
pseudomomentum deviations as 

j = j<^ + j^ + j^,,,^ (A.3) 

where the first term, J^, of the current normal-ordered expansion (|A.3|) can after 
some algebra be written as 



J' 



j=±i 



sin(Q) + 



2vrpc,o(Q) 

7>o j=±i y 1 - M2[jrc,7 - 



7 V c,7 ) 
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where the functions 27rpc,o(^) and 2n pa^^{r) were defined in Ref. and 



+ }_^ReAt dqSNcM^ '\ = — f — L^'^iq) ■ (A.5) 



7>0 



The function jC^'^{q) is defined as 

C''Hq) = L''%)-W\q)^^. (A.6) 

In order to obtain the integral equations for L'^'^lq) and C^''^{q) (with C = 
-^c.G, -^£,7, and Cs^-y), we start from the continuum hmit of Eqs. ( |6.9|) , ( |6.10|) , and 
( |6.11| ) which reads 



V>o ^-9c,y V Y 7 



2i?e sm-^[{RcM - i-f)u] =q + 7(0^ + <f)i)/Na - 



— r dq'N,{q') 2 tan 

ZTX J-Pr 



7'>0 



and 



- 'i'(XT('(')ei+i,T+i(fl.,7fe)-fl.,T('('))- (A.9) 

, /TT J-qs,-i 



7' 



It is convenient to write the function 0'^L(a;), defined by Eq. ( |A.9| ) of Chapter 



as follows 
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We emphasize that comparision term by term of Eq. (|A.9|) of Chapter ^ with 
expression ( |A.10| ) above fully defines the coefficients b]''^ and the corresponding set 



of integer numbers /. 

Following equation ( |6.18|) , we have that 0f = 0^ for a charge-probe current and 



(pl = —(pi for a spin probe. With the above equations written in terms of 0| and (pi, 
Eq. ( [A.4| ) provides both the charge and spin currents. In what follows, we introduce 
in the functions L'^{q) the index ( = p,az to label the equations for either the charge 
or the spin current, respectively. We start by expanding Eqs. ( |A.7| ), (|A.8| ), and ( |A.9[ ) 
up to first order in 0. This procedure reveals that the functions L'^''^{q) obey the 
following integral equations 



A^,,y(gO dRs,y{q') ^^^^ , 



+ 



r sin(X(g))/n-fl^y (gp p srt' 



y+1 



^c^Q) - ^c,y + J_^^ dq ^ ^ ^ 3in(i.(,0)/.-i^.,,(,) p cos(K(g jjL,,o(? j 

and 



s,7 



+ 



^(7 + VjTi 

7 ' 



dRs^yiq') 



s,7 



(A.13) 



where the coupling constants C^,^ are defined by Eqs. ( p. 25] ) and ( p.26[ ). We 
again write the distributions functions Na^.y{q) of Eqs. ( |A.11|) , ( [A.12|) , and (|A.13| ) 
as N^^^iq) + SNa^-yiq). This allows us to obtain integral equations for L^''^'^{q) and 
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£^''^{q) (we remark that the functions £^''^{q) are the same both for ( = p, o"^). It is 
then straighforward to find the integral equations obeyed by L^''l''^{q) and show that 
L^''^'^{q) can be simply expressed in terms of linear combinations of phase shifts. 
The final result is 

= Ci, + E E J^(A^a',y)CS,,,$„,^;a.,y(g, . (A.14) 

a', 7' i=±l 

The integral equations obeyed by C^''^{q) arc related to the integral equations obeyed 
by £^''^(r), where r equals sin(_ft'('^)(g))/-u, Rf^^{q), and Rf}y{q) for C = £c,0) -^0,7, 
and Cs,'y, respectively. The functions C^''^{r) obey the following integral equations 

and 



1 pre f^''^(r') 

£},<p{r] = C^''t>'^(r) = / dr' "'^^ ' 

J~'s,^r) L.,^^ ^rj ^ ar ^ ^ ^^^^ 

where the free terms £c'o'°(r), >Cj'^'°(r), and >C^'^'°(r) are, respectively, given by 



E-7 



1 fJe.y 



dg'(57V^,y(g')- 



l + [ 



c,7 12 



(i + l^P) 



+ E^(^c,7')-;- E 



1 V- i^c,7'(i?Fc,7') L^'UjqFc. 



+ E^(^^^ 



s,7 



E 
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1 



■nwy J-qc 



dq'6N,{q')- 



1 + 



|- sin(i^(°)(g'))/-u-r i2 



cos(A'<«)(,'))-'''^"°'"'> 



7' ' 



l + [ 



1 V- JLIoUqfc) lI'^U^Ifc) 
^UT^^ii, 2vrp,,o(g) (1 + 



' J=±l 



27rpc,7'(^c,7' 



(l + [^]^) 



dq' 



(A.19) 



and 



™(7 + 1) i-,. ^ ^ 1 + [ sin(i^W(g'))/n-n 2 



7' i '^s,7' ]^ 



sin; 

7+T 
7+l,7'+l rO, 



dq' 



hr'''^'L';Uq')dR':Uq') 



+ 



1 



(7 + l)«vr ^.f^^ 27rp,,o(Q) (1 + [^l^) 



EW,y)E:j7E- 



7+l,7' + l r 1 



2vrp, 



s,7' V s,7'> 



(1 



(A.20) 



Introducing the functions C^'^{q) obtained from Eqs. ( A.15 ), ( A.16|) , and (|A.17|) 
in Eq. ( A.4 ) and keeping terms only up to second order in the density of heavy 
pseudoparticles, we obtain Eq. ( |6.19|) with ji^^iq) given by 



£,7(5) = ^a,7(g)^°'i'^(9)+ E E j0iNc',y)Va',^'LyY{jqFa',y)^a',y;a,^{jqFa',y,q) ■ 

a' ,7' j=±l 

(A.21) 

Inserting Eq. {J^J§ in Eq. {^^^) we obtain Eq. ( ^ ). 
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Appendix B 



Static masses for the heavy 
pseudoparticles 



The a, 7 pseudoparticle static mass, m* , is defined in Ref. as 



2t dria^-y{r) / dr 



m 



0,7 



(2vrp„,^(r))2 



2tr]a,y{r) (27r dpa^yir) / dr) 



fB.r 



=^0 (27rp«,^(r))3 
where the functions 2trja^^{r) and 27rpQ .^(r) are defined in Ref. and r° stands for 
^arMPcin) which represents Q, r^^, and r^^-y. 

After some straightforward algebra, the general expressions ( [B.l| ) lead to the 
following simple expressions for ^ 



c,7 



and 



-4tuV(l + ^272)3/2 + A^^^ 



A*? _ - A''' 



7>0: 



s,7 



In Eqs. (p.2| ) and ( |B.3|) the functions AJ^ ^ and A^ ^ read 



7 > 0. 



QF 



and 



A^ 



qpa.o dq 



1 



,0 TTX 1 + o(?)/a;]^ 



with 



fB.2l 



(B.3) 



(B.4) 



(B.5) 
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In the limit of fully polarized ferromagnetism, these expressions lead to the 
following closed-form expressions for the static masses 



1 tn I TT + 2[?7i^^] U7sin(2n7r) 



+ sin (mv) 



and 

1 tn I 1 ui^j + 1) sin(2n7r) 



:b.7) 



:b.8) 



rn^-y [^2,7+i] + [m(7 + 1)]2 2[r/2,7+i] [u(7 + 1)]^ + sin2(n7r)^ 

where 

7]!,^ = tan~^[cot(n7r)^===] , (B.9) 
V 1 + M a;^ 

and ri2^x = 7r/2 — r/i -j,. We remark that the static masses of the c, 7 pseudoparticles 
are, in general, negative. The static masses of the a, pseudoparticles have been 
studied in Ref. [HS . 



Appendix C 



Charge and spin stiffnesses at zero 
temperature 



In this Appendix we show that the direct use of Eq. ( |6.35| ) leads to the stiffness 
expressions (135) - (137) of Ref. 

The calculation of the charge and spin stiffnesses (|6.35|) requires the expansion of 
Eq. ( |6.12| ) and of Eqs. ( |A.7| ), ( |A.8| ), and ( |A.9| ) up to second order in 0. As in the case 
of the charge and spin current, both he charge and spin stiffnesses can be computed 
from Eq. ( |fci.35D , and we obtain one or the other depending on the coupling constants 
we choose in Eqs. ( |A.7| ), ( |A.8| ), and ( [A.9| ). Expanding the ground-state energy up 
to second order in ch, we obtain 



d{<PlNa) 



1 flp' 



(/>=0 



27r 



dq [2tK°''^<^(g) sin(ir(o)(g)) + 2t[ir°'<^(g)]2cos(ir(°)(9))l 



(C.l) 

where the function K^''^'^{q) is the second derivative of the rapidity function defined 
by Eq. (pTZl ) in order to at = 0. The functions K^^'l'{q) and K^^'l"t'{q) can 

be written as 



dq 



and 



(C.2) 



K'^^^^q) = ^ 



d (dK('Kq)no,H^,A , ^dK('Hq)^o,H 



■[^c;oW +2- 



c,0,* ) 



dq \ dq ^ ^'^'^ J ' dq 
respectively. The use of Eqs. (|C.2|) and ( |C.3| ) in Eq. (|C.1|) leads then to 



(C.3) 
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^ ^r.,2tsm(ir(0)(,))2^Lat(,) 

27r 27rp,,o(Q) ' ^ ' ^ 



where the function L^ct^JlFcj is defined in Appendix ^ The function L°'ot(?) obeys 
the following integral equation 



27r ^tii 27rp,,o(r.,o)(l + [sin(ir(o)(g))/n - jr^flf) 

7ri-,,„o ^ dq' l + [sm{Ki^){q))/u-R%{q')]^' 

which was obtained by performing the type of expansions developed in Appendix 
0. Moreover, is given by 



1_ j[cosiQ)Ll^{jqFcW 
uvr^.f^, 27rp,,o(Q)(l + [i??^(g)-J 
1 „ j[L''tijqFs,o)r 







E 



1 r^^^o ^^^dR^^liq') L°;Jt(g') 



^ ^ ^ , dK(^\q') cos(irW(gO)L°ff(gO 

TTwi-g^, ^ rfg' l + (sin(ir(0)(g'))/M-i?f,i(g))2 ■ 



Introducing Eqs. (|C.5|) and (|C.6|) in Eq. ( p.4|) we obtain, after some algebra, the 



following expression for the (charge and spin) stiffness Z)^ 

AnD^ = Y: v,,o[L'SHjQfc,o)? + E ^.,o[<o^''^(jgFs,o)]' , (C.7) 
j=±i j=±i 

where the functions L^''^'^{iqFafl) are defined by Eq. ( |A.14| ). After some simple 
algebra, expression ( |(J.7D can be shown to be the same as expressions (135) - (137) 
ofRef. 



Appendix D 

The zero magnetic-field solution 



For the case of zero magnetic field it is possible to cast the equations for the energy 
bands, phase shifts, and rapidities in a simpler form. After some algebra, the €^^{q) 
band (with 7 = 0, 1, 2, ... , 00) and the e°^(g) band with (7 = 1, 2, ... , 00) can at 
zero magnetic field be rewritten as 



and 



2t I duj , ' . Ti(^) 

a;cosh(u;j 



e° (g) = ReU^l-u^iR^^liq)+t^y- 



4t / duj cos(cui?(°)(g))Ti(cu) , 

JO UJ " 



where Ti{uj) obeys the integral equation 



(D.i: 



(D.2) 



/oo 
rfcu'Ti(cu')r(w', 
-00 

Here the free term and the kernel read 



(D.3) 



2tt J-q 



dk sm{k) sm{uj sin( A;) / u) 



and 



(D.4) 



sm{{uj - uj')rc,o) 



:d.5) 



7t{uj - io'){l + e^\^'\) ' 
respectively. The kernel ( p.5| ) was already obtained in Ref. (see Eq. (A5) of 



that reference). Equation (|D.1|) , together with the fact that in the limit of zero 
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magnetic field the width of the s, 7 > momentum pseudo-Brillouin zone vanishes 
[see Eq. ( |6.36|) ], shows that the s, 7 bands collapse for 7 > and all values of U and 



n to the point zero. 

In this limit it is also possible to cast the integral equations for the phase shifts, 
whose expressions are given in Ref. in the following alternative form 

$c,0;c,o(r, r') = -B{r - r') + f"' c/r"i>,,o;c,o(r", r')A{r - r") , (D.6) 



$c,0;c,y r') = -- t^n\'-^) + dr"<l>,,o;c,y (r", r')A{r - r") , (D.7) 



^c,Q;sa'i.f,r') = -5o,y7^tan ^[sinh(7r/2(r -r'))] + 

Zn 



+ r -^^c,0;s,Y{r", r')Air - r") , (D.8) 

J-Tcfi 7' + 1 



$c,7;c,o(r,r ) = -tan ( )-/ ,r-r"N9 ' (^-9) 



TT 7 J-r.fi 7r7 1 + (^) 

$ ,(r r') = ,(r - r') - ^ i£:2i£>:!lll (d 10) 

^c,7;c,yl^,n 27r^'^ ^ ^ i-r.o 7r7 1 + ' ^ ' 

5's,-y;c,o(r,r') = 0, (D.12) 



$ ,fr r') -5, - r dx" fD 13) 



sin(u;[r — r']) 



TT Jo u;(l + e^"^) 



[e-^'" + (l-5o,y)e 



-(7'-2). 



1 r<=>0 dr" $c,0;s,7'(^")^') 

^ 4 y-r,,o 7' + 1 cosh(7r/2(r-r")) 



(D.14) 



and 
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*s,7>0;s,7'('^,'^') = 77-07+1,7'+! l'" ^ ^ 

The functions yl(r) and B{r) are defined as 



TT JO 1 + e^''^i 



and 



B(r) = lrd. (D.17) 

respectively. 

At = 1 we have that Ti(a;) = T^{x) = Ji{x/u), where Ji{x/u) is the Bessel 
function of order one, and the bands ( p.l|) and (p.2| ) are obtained in closed form. 
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Chapter 7 
Conclusions 



In this Thesis we have focused our attention on the one- dimensional Hubbard model. 
As it became clear from Chapters |l] and the Hubbard Hamiltonian, as introduced 
by Hubbard, Gutzwiller, and Kanamori is the simplest model for the 

description of interacting electrons in a lattice. 

Introduced firstly to study itinerant ferromagnetism in transition metals, the 
Hubbard model revealed later to have an antiferromagnetic ground state [^. With 
the discovery of a new class of superconducting materials, by Berdnoz and Miiller in 
1986 and the proposal, by Anderson that the 2D Hubbard model would retain 
the essential physics of these new materials, the model gained renewed interest. It 
soon became clear that, for the values of the physical parameters imposed by the 
experimental results, the usual theoretical tools would not suffice to deal with this 
model (and its "relatives") 0] and that, accordingly, the search for non-perturbative 
methods should be pursued. The claim by Anderson that the 2D Hubbard model 
should present many similarities with its one-dimensional counterpart , for which 
an exact solution is available, together with better crystals of quasi one- dimensional 
materials, and more refined experimental methods (e.g. angle resolved photoemis- 
sion measurements) P] , renewed the interest in the one- dimensional Hubbard model. 

At present, the theoretical goal associated with the study of the Hubbard model 
is to investigate the role of electron-electron interactions in relation with the ther- 
modjTiamic and transport properties of a non-perturbative electron liquid. Based 
on the exact Bethe-ansatz solution of the Hubbard chain and on its symmetries, 
we presented an algebraic representation for all the eigenstates of the model (Chap- 
ter 0). Our representation is much simpler than the starting Bethe-ansatz wave 
function. The latter is a complicated sum of permutations written in terms of 
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the original electronic operators, while the former representation is a simple Slater 
determinant of pseudoholes and heavy pseudoparticles (Chapter H). These latter 
operators are the quantum operators that diagonalize the Hubbard Hamiltoninan. 
In this representation it is then simpler to select and characterize the relevant eigen- 
states which contribute to a given physical property of the model. For example, 
the low-temperature thermodynamics for electronic fillings less than one and finite 
magnetization is controlled by the eigenstates with no heavy-pseudoparticles oc- 
cupancy; on the other hand, the optical conductivity is controlled, at half filling, 
by eigenstates with one c, 1 heavy pseudoparticle and two c, pseudoholes. This 
kind of characterization is by no means possible for the usual Bethe-ansatz wave 
function, and it is essencial for a simple calculation of correlation functions using 
conformal-field theory methods ^ 10, |ll|, [T^ . 

The pseudoparticle perturbation theory (PPT) described in Chapter |, combined 
with the results of Chapter and with new results on conformal theory methods 



111 , 0, allow the calculation of correlation functions at low {u — uq) frequency 
values (note that Uq can be very large). This possibility goes much beyond both 
bosonization studies |T^, and the spinon and holon conformal-field theory of Framh 
and Korepin @, [l^l- The possibility opened up by our pseudoparticle theory, namely 
that of giving finite energy results for the spectral functions of the model, allows 
the comparison of finite energy spectral functions in real materials with the results 
obtained for the Hubbard model. This is a topic that deserves further research. 

The use of the algebraic representation introduced in Chapter ^ permits the 
computation of the energy eigenvalues for all the eigenstates of the model. This can 
be accomplished, at least in principle, by means of the PPT developed in Chapter 
^. This PPT is an expansion of the Hamiltoninan in the density of excited pseu- 
doparticles, relatively to a suitable reference state. The most convenient choice for 
the reference state is, for most of the physical properties, the ground state or the 
generalized ground-state, both introduced in Chapter ||. 

The normal-ordered PPT Hamiltonian describes the physics associated with the 
Hilbert subspace spanned by the reference state and all Hamiltonian eigenstates 
differing from it by the occupancies of a small density of pseudoparticles. Com- 
bined with a generalized conformal-field theory, the PPT provides the leading-order 
expressions of correlation functions for low-energies {uj — ujq)- These expressions 
refer either to the low-energy functions (tuo = 0) [10|, Q or to energy values 
closed to the energy gaps, ujq, of Hamiltonian eigenstates with heavy-pseudoparticle 
occupancy [Tl|, |T^. 
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Very important for the description of excitations that change the number of 
electrons A^^ and/or (change of canonical ensemble) and/or the number of pseu- 
doparticles (change of sub-canonical ensemble) is the study of the collective pseudo 
momentum shifts of value ivr/iVa. These momentum shifts are generated by the 
topological momentum shift operators introduced in Chapter |[ These are impor- 
tant in two ways: (i) they contribute to the momentum of state transitions con- 
necting different canonical ensembles and sub-canonical ensembles [p!5| , |I6|; (ii) in 
transitions between ground states, the final pseudoparticle momenta are such that 
the final ground state has the lowest possible energy [l^ . The importance of these 
types of excitations had not been recognized before in the literature in its full rele- 
vance. 

The problem of charge and spin transport in the one-dimensional many-body 
problem, in general, and in one-dimensional integrable models, in particular, is a 
difficult and subtle issue |T^, WU[ EO, pTI, By using the Bethe-ansatz equations 



with twisted boundary conditions, we have obtained a general formula for the mean 
value of the charge and spin current operators (Chapter However, the functional 
character of this expression does not seem amenable to practical calculations. Com- 
bining this general functional with the PPT, and performing equivalent calculations 
to those presented in Chapter ^, we were able to compute the mean value of the 
charge and spin currents for a given subset of eigenstates of the model. The ob- 
tained normal-ordered functional provides the current for states that differ from the 
reference state by a small density of excited pseudoparticles. Our current functional 
provided the values for the transport mass of the charge and spin carriers. These 
masses provide important information on the role of electronic correlations in the 
transport properties. Our functional representation can also provide information on 
the transport of charge and spin at finite temperatures. Althougt the PPT does not 
cover the full range of temperatures T from to oo, it provides some information 
in some regimes. For example, away from half filling, the regime where PPT can be 
used is T ^ vfcaQfca-, while at half filling it is T <^ A, where A is the smallest of 
the energy gaps c^o derived in Chapter ^. Studies on finite-temperature transport 
by use of the PPT is a subject of current research. 

The study of the Bethe-ansatz equations, together with the PPT, gave the 
coupling constants of the pseudoparticles to the external charge and spin probes. 
These coupling constants provide selection rules for the ground state transitions. 
Conformal-field theory provides asymptotic expressions for correlation functions, 
including their x, t (space and time) and uj (momentum and frequency) depen- 
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dences. However, it does not provide the multiplicative constants determined by the 
weights of the ground-state transitions, which control these expressions and depen- 
dences. Often, due to symmetry and selection rules, these weights and corresponding 
constants vanish. Therefore, the above couplings provide relevant information on 
such rules, which goes beyond conformal-field theory. 

Both the coupling constants and the mean value of the charge and spin current 
operators can also be derived from a semi-classical transport description of the 
quantum liquid. Comparing this approach to the Bethe-ansatz exact solution (PPT 
theory), one confirms that the kinetic approach gives exact results for small {uj — ujq), 
while for a Fermi liquid the same kinetic approach gives exact results for small u 
only. This difference can be traced back to the existence of only forward scattering 
among the pseudoparticles at all energy scales, and it is peculiar of integrable one- 
dimensional models. 

The formalism presented here is applicable to other one-dimensional integrable 
systems Examples of these are the super-symmetric t—J and the lattice spinless 
fermion models. Both these models present interesting transport properties, since 
the charge current operator does not commute with the Hamiltonian. Besides the 
charge stiffness, there is a finite-energy absorption that can be studied by combinig 
the methods presented in Chapters |^, ^, ^, and |^ with the thechniques of conformal- 
field theory p], |10|, |Tl], |12|. Also the study of finite-temperature transport can be 
fulfilled by developing the apropriate PPT for these models. Following the results 
presented here, calculations for the lattice spinless fermion model are currently being 
developed [|4 . 
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